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ABSTRACT 


A  submerged  gold  wire  was  electrically  exploded  to  simulate  a  very 
deep  underwater  nuclear  detonation.  These  small-scale  tests  were  con¬ 
ducted  in  a  test  tank  in  which  the  air  pressure  may  be  varied  to  obtain 
various  strengths  of  buoyant  bubble  migration. 

The  mechanics  of  explosion  bubble  formation,  pulsation,  and  migra¬ 
tion  are  discussed,  and  relevant  relationships  are  derived  far  scaling 
large-yield  bubble  behavior  in  the  laboratory  by  variation  of  environ¬ 
mental  parameters.  The  mechanisms  by  whicn  the  explosion  products, 
initially  at  the  bubble  center  during  the  first  cycle,  are  redistributed 
by  bubble  bottom  collapse,  compression,  re -expansion,  and  migration, 
are  investigated. 


The  experimental  technique  involved  physically  sampling  the  bubble 
and  its  environment  at  the  second  bubble  maximum  to  define  the  extent 
and  nature  of  explosion  product  re-distribution.  Three  different  shot 
conditions  were  investigated,  corresponding  to  a  "rest-point"  condi¬ 
tion,  a  condition  characterized  by  slight  upward  migration,  and  one 
characterized  by  strong  upward  migration. 


In  all  cases,  it  was  found  that  most  (77  to  9k%)  of  the  explosion 
products  are  lost  to  the  environment  by  the  time  of  the  second  bubble 
maximum.  The  probable  mechanisms  for  this  loss  are  discussed,  and  rules 
are  derived  for  application  of  these  results  to  full-scale  under¬ 
water  bursts. 
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SUMMARY 


The  Problem 

When  a  nuclear  explosion  occurs  underwater,  some  of  the  radio¬ 
active  fission  products  are  lost  to  the  liquid  environment  as  the  steam 
bubble  pulsates  and  migrates  upward.  The  remainder  are  ejected  into  the 
above-surface  environment.  In  order  to  quantify  the  extent  of  thiB  cub- 
surface  loss,  it  is  first  necessary  to  identify  and  evaluate  in  detail 
the  mechanisms  by  which  explosion  product  material  is  transferred  from 
its  initial  location  at  the  center  of  the  bubble  to  the  environmental 
water  by  the  hydrodynamic  processes  occurring  during  bubble  collapse, 
migration,  and  re-expansion.  The  objective  of  this  study  is  to  evaluate 
these  phenomena  as  a  function  of  the  strength  of  bubble  migration. 

o 

Findings 

In  the  case  of  no  migration,  only  about  185 i  of  the  explosion  pro¬ 
ducts  were  found  within  the  bubble  at  its  second  maximum.  If  weak  up¬ 
ward  migration  occurs,  the  retained  fraction  drops  to  about  6$,  but,  if 
migration  strength  is  further  increased,  the  underwater  Jet  and  de¬ 
creased  internal  energy  of  the  bubble  at  the  minimum  cause  less  material 
to  be  lost  -  the  fraction  remaining  in  the  bubble  at  the  second  maximum 
for  the  "strong  migration"  case  is  about  23$  of  the  total. 
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CHAPTER  1 


INTRODUCTION 

With  the  development  of  nuclear  weapons  for  use  in  antisubmarine 
warfare,  contamination  of  the  above-surface  environment  by  the  radio¬ 
active  fission  debris  from  the  weapon  has  become  a  problem  of  some  con¬ 
cern.  Even  if  the  burst  is  quite  deep,  seme  fraction  Gf  this  residue 
will  be  carried  to  the  surface  by  the  migrating  explosion  bubble,  and 
will  become  involved  in  the  plumes  and  the  radioactive  base  surge  which 
may  represent  a  serious  hazard  to  surface  vessels. 

It  has  always  been  assumed  that,  if  the  depth  of  burst  is  suf¬ 
ficiently  deep  so  that  the  explosion  bubblr-  will  pulsate  several  times 
before  reaching  the  water  surface  (the  "very  deep"  burst  classification),^ 
a  substantial  fraction  of  the  fission  products  will  be  transferred  from 
the  bubble  to  the  water  environment,  thus  reducing  the  above-surface 
hazard.  Quantitative  information  concerning  this  loss,  .  ever,  is 
not  available.  There  has  been  only  one  full-scale  nuclear  detonation 
falling  into  the  "very  deep"  category  (Operation  Wigwam),  and  the  data 
obtained  concerning  post -shot  fission  product  distribution  was  too 
fragmentary  and  incomplete  to  evaluate  the  relevant  loss  mechanisms. 

It  is  intrinsically  quite  difficult  to  estimate,  for  example,  the  frac¬ 
tion  of  the  total  activity  present  in  the  plumes,  both  because  radiation 
field  data  is  rather  poor  and  the  water  density  structure  of  the  plume 
is  not  known.  An  attempt  was  made  at  Wigwam  to  determine  the  total 
amount  of  activity  present  in  the  subsurface,  radioactive  "pool",  but 
measurements  were  not  possible  at  early  times,  and  by  the  time  radio¬ 
activity  was  encountered,  oceanographic  processes  had  so  distorted  and 

p  q  4 

diffused  the  debris  that  no  reliable  result  was  obtainable.  In 

any  case,  the  vertical  distribution  of  fission  products  in  the  water  to 
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be  found  several  days  after  burst  is  not  necessarily  characteristic  of 
that  prevailing  at  early  times. 

Furthermore,  simulation  of  the  radiological  transfer  processes 
for  deep  and  very  deep  bursts  by  means  of  traced,  relatively  low  yield 
high-explosive  field  tests  is  not  desirable,  for  two  reasons.  First, 
the  strength  of  bubble  migration  is  yield  dependent,  and  consequently 
approximations  must  be  made.  Second,  the  explosion  bubble  atmosphere 
in  the  high-explosive  case  consists  of  non-condensable  gaseous  reaction 
products,  rather  than  steam  as  in  the  nuclear  case.  Developmental  work 
is  now  in  progress  on  high-explosive  charges  producing  steam  as  a 
principal  reaction  product,  but  such  explosives  are  not  as  yet  available. 
The  importance  of  a  condensable  bubble  atmosphere  to  explosion  product 
transfer  processes  will  be  discussed  at  length  in  later  chapters . 

The  conceptually  simple  approach  of  firing  a  large  explosion  and 
then  taking  measurements  after  violent  motion  has  subsided  may  not 
prove  particularly  fruitful  in  any  case.  In  order  to  generalize  results, 
more  than  an  a  posteriori  distribution  for  some  particular  charge  fired 
at  some  particular  depth  at  some  particular  time  after  burst  is  requir¬ 
ed.  That  is,  a  prerequisite  to  any  method  of  prediction  of  the  distri¬ 
bution  of  activity  is  an  understanding  and  evaluation  of  the  processes 
by  which  material  is  lost  to  the  environment  by  the  bubble  as  it  pul¬ 
sates  and  migrates  to  the  surface,  as  functions  of  yield  and  shot  depth. 

For  this  investigation,  the  experimental  work  was  carried  out  in 

a  vacuum  tank  using  an  electrically  exploded  gold  wire,  which  is  a  good 

simulant  to  a  nuclear  explosion  since  it  generates  a  steam  bubble  in 

which  the  explosion  products  are  initially  distributed  in  the  appro- 
5 

priate  manner.  In  order  to  simulate  large-yield  bubble  behavior 
(migration  strength)  the  air  pressure  over  the  water  was  varied.  Three 
shot  conditions  were  investigated;  one  for  which  there  was  no  migration. 


2 


one  for  which  migration  was  unward  but  relatively  weak,  and  one  for 
which  upward  bubble  migration  was  a  dominant  feature  of  the  motion. 

The  initial  distribution  at  the  first  bubble  maximum  lias  previously 
been  documented.^  In  the  present  work,  in  order  to  identify  and  eval¬ 
uate  the  hydrodynamic  mechanisms  by  which  the  explocion  products  are 
re-distributed  in,  and  ejected  from,  the  bubble  during  bubble  collapse, 
translation  and  re-expansion,  the  bubble  and  the  fluid  environment  from 
each  of  these  three  conditions  were  physically  sampled  at  the  time  of 
the  second  maximum.  The  samples  were  then  analyzed  for  explosion  pro¬ 
duct  (gold)  content  by  neutron  activation  techniques.  From  these  dis¬ 
tributions  of  explosion  products,  processes  are  postulated  based  on 
knowledge  concerning  the  flow,  both  inside  and  outside  the  bubble,  near 
the  time  of  recompression. 

A  prerequisite  to  any  understanding  of  the  mechanisms  of  loss  and 
re-distribution  of  explosion  products  is  familiarity  with  the  bubble  phe¬ 
nomenon  and  the  type  of  fluid  motion  to  be  expected  after  an  underwater 
explosion.  Consequently,  Chapter  2  deals  with  the  present  state  of  the 
art  concerning  bubble  mechanics  and  scaling  of  bubble  behavior.  Later 
chapters  will  discuss  such  matters  as  the  experiment  itself,  data  re¬ 
duction,  and  quantitative  results  concerning  the  mechanisms  of  explo¬ 
sion  product  loss. 
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CHAPTER  2 


PHENOMENOLOGY  AND  SCALING  OF  UNDERWATER  EXPLOSION  BUBBLES 

FORMATION  OF  THE  BUBBLE 

The  sequence  of  events  following  a  deep  underwater  explosion  depends 
to  some  extent  on  the  nature  of  the  explosive.  Generally  speaking,  we 
may  divide  explosives  into  two  major  classes.  First,  we  have  conven¬ 
tional  high-explosives  which  liberate  energy  by  high-speed  chemical 
reactions  and  have  energy  densities  less  than  3 >000  calories/cnr  or  so. 
The  "point-source"  explosives,  on  the  other  hand,  are  characterized  by 
energy  densities  at  least  an  order  of  magnitude  higher  than  those  of 
most  high-explosives,  and  therefore  have  non-chemical  energy  sources; 
these  include  high-energy  electrical  discharges,  electrically  exploded 
conductors,  and  nuclear  devices. 

Considering  first  the  high-explosive  case,  when  the  detonetion  wave 
reaches  the  surface  of  the  charge,  it  proceeds  into  the  water,  leaving 
behind  it  a  sphere  of  gaseous  reaction  products  at  high  temperature  and 
pressure.  After  this  time,  the  interaction  between  the  shockwave  and 
these  residual  gases  may  be  considered  negligible.  It  may  be  reasonably 
assumed  that,  at  this  stage,  this  "initial  bubble”  is  relatively  homog¬ 
eneous  as  regards  temperature  and  density.  These  gases,  as  a  conse¬ 
quence  of  their  high  pressure,  do  work  on  the  environment  and  expand 
adiabatically  in  a  spherically  symmetric  fashion. 

In  the  case  of  the  point-source  detonation  with  its  much  higher 
energy  density,  the  materials  constituting  the  device  are  immediately 
raised  to  extremely  high  temperature,  and  are,  in  general,  ionized. 

The  emission  of  high  energy  photons  and  the  radiation  of  heat  from  the 
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device  causes  vaporization  of  the  water  in  the  vicinity.  The  resulting 
high  pressure  system  expands  as  a  shock  front.  As  this  shockwave 
passes  through  the  adjacent  water,  the  water  is  first  raised  to  high 
pressure  and  temperature  above  the  critical  point,  and  then,  as  the 
shock  front  passes  on,  the  pressure  drops  and  the  water  goes  into  the 
vapor  phase.  At  fairly  early  times,  the  shock  front  surface  energy 
density  drops  below  that  value  required  to  further  vaporize  water,  and 
separation  between  the  bubble  front  and  the  shock  front  occurs.  As 
can  be  seen,  the  two  mechanisms  of  steam  generation  (initial  radiation 
and  shock  energy  dissipation)  both  tend  to  deposit  less  energy  per  unit 
mass  of  water  with  increasing  distance  from  the  device.  Consequently, 
in  the  point-source  case,  the  gas  temperature  is  highest  and  the  density 
is  lowest  in  the  center  of  the  initial  bubble.'1' 

The  volume  of  the  bubble  at  the  time  of  separation  is  about  the 
same  as  that  of  a  high  explosive  charge  of  the  same  yield.  Hence, 
after  this  time,  the  initial  bubble  motion  should  be  similar  for  both 
cases.  However,  the  fact  that  the  high-explosive  bubble  atmosphere  is 
identical  with  the  explosion  products,  whereas  the  debris  remaining  of 
the  point-source  device  is  located  at  the  center  of  a  matrix  of  condens¬ 
able  steam,  causes  appreciable  differences  in  behavior  at  later  times, 
.■tore  will  be  said  about  this  later.  The  best  qualitative  distinction 
between  a  point -source  explosion  and  the  usual  high-explosive  detona¬ 
tion  may  therefore  be  made  in  terms  of  the  composition  of  the  bubble 
atmosphere,  which  is  in  turn  related  to  the  energy  density  of  the  explo¬ 
sive.  That  is,  high-explosive  energy  densities  are  typically  in  the 
range  1000-3000  calories/ cur*;  it  is  probably  safe  to  assume  that  be¬ 
havior  will  be  characteristically  "point-source"  (that  is, the  bubble 
atmosphere  will  be  primarily  steam)  if  the  explosive  energy  is  greater 
than  10,000  cal/cm^. 
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THE  NON-MIGRATING  CASE 


For  the  present;,  in  order  to  simplify  the  initial  development, 
the  effects  of  gravity  will  be  neglected,  and,  furthermore,  it  will  be 
assumed  that  both  the  sea  surface  and  bottom  are  far  away  from  the 
explosion  point.  Under  these  circumstances,  the  bubble  remains  spher¬ 
ical  and  the  problem  is  one -dimensional.  The  more  generei  case  will 
be  taken  up  later  in  this  chapter. 

After  the  initial  bubble  has  been  formed,  and  the  shock  wave  has 
separated  from  the  bubble  front,  the  motion  for  high-exploBive  and  point- 
source  bubbles  is  almost  identical  until  much  later  times.  By  virtue 
of  the  high  internal  pressure,  the  system  expands  radially,  initially 
with  very  high  velocity.  As  the  bubble  does  work  on  the  environment, 
generating  kinetic  energy  in  the  adjacent  water,  the  internal  energy  of 
the  adiabatically  expanding  gas  drops  rapidly.  At  a  fairly  early  time, 
the  internal  pressure  becomes  equal  to  the  external  (hydrostatic) 
pressure.^  At  this  time,  most  of  the  total  energy  resides  as  kinetic 
energy  in  the  water  environment.  After  this  time,  the  environment 
begins  to  do  work  on  the  bubble  (the  net  force  is  now  inward),  eventu¬ 
ally  bringing  the  expansion  to  a  stop,  much  later.  Then,  the  system 
begins  to  contract,  and  the  process  is  reversed.  The  bubble  internal 
pressure  increases  to  some  peak  value,  and  the  system  re-expands. 

These  cycles  continue,  with  decreasing  amplitude  and  period,  until  all 
the  available  energy  has  been  lost. 

The  maximum  radius  and  the  period  of  the  bubble  for  the  first 
cycle  are  related  to  the  total  energy  released  and  to  the  local  hydro¬ 
static  pressure.  In  order  to  avoid  ambiguities  concerning  explosive 
type,  which  lead  to  correction  factors  which  must  be  carried  throughout 
the  development,  the  basic  unit  of  explosive  bubble  energy  will  be 
taken  as  the  ''pound -TNT-equivalent -bubble".  This  transformation  allows 
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the  rest  of  the  development  to  be  done  in  terms  of  dimensionless 
quantities.  If  the  weight  of  the  charge  (which  is  proportional  to  the 
total  energy  available,  for  a  given  explosive)  is  denoted  by  W*,  then 
we  define  the  "TNT  equivalent  bubble  yield"  as  follows: 

W  =  2W*  'e/eTJ  ,1  -  fs 

J  L  V 

wuere 

W  =  TNT  equivalent  bubble  yield  -  expressed  in  pounds 
e  =  cp/pe  f°r  the  explosive  of  interest 
CTNT  =  q)/pg  for  TNT 

cp  =  energy  density  of  the  explosive 
Pe  =  density  of  the  explosive 
^  =  total  energy  released 

Es  =  energy  initially  radiated  in  the  shockwave 

It  can  be  seen  that  the  term  (l  -  (Eg/E^))  simply  represents  the 
fraction  of  the  total  energy  released  which  remains  available  for 
bubble  motion.  This  quantity  does  not  change  with  yield,  and  is  char¬ 
acteristic  of  the  explosive  used  -  the  appropriate  value  for  THT  is 
about  0.5.  If  the  energy  density  of  the  explosive  is  high  enough  (as 
for  a  point-source  explosion),  losses  may  also  occur  through  electro¬ 
magnetic  radiation,  in  which  case  the  above  term  should  be  replaced  by 
(l-((Eg  +  Ep)/!^,)).  The  value  of  may  be  taken  as  1000  calories/gram. 
We  may  then  simply  say  that 

W  =  kW’  (2.1) 

where  k  is  characteristic  cf  the  explosive.  Values  of  k  for  several 
common  explosives  are  listed  in  Table  2.1. 
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TABLE  2-1 


Values  cjf  the  Bubble  Yield  Conversion  Factor  for 
Various  Explosives 
(Compiled  from  various  sources) 
k  =  W/W' 


s 

k 

Amatol  (cast) 

1.1 

HBX-1 

1.5 

Lead  Azide 

0.36 

Minol-2  (cast) 

1.4 

Minol-2  (pressed) 

1.7 

Pentolite 

1.0 

PETN 

1.6 

HDX  (loose  mix) 

1.2 

BDX  (pressed) 

1.0 

Tetryl  (loose) 

1.1 

Tetryl  (pressed) 

1.0 

TNT 

1.0  (by  definition) 

Torpex 

1.5 

Exploding  vire) 

Spark  Gap  ) 

0.6  -  1.0  (depends  on  circuit  parameters, 
etc.) 
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The  hydrostatic  pressure  at  the  shot  point  is  denoted  by  Z  and  is 
expressed  in  feet  of  water.  Hence,  if  the  air  over  the  water  is  at  one 
atmosphere  pressure, 

Z  =  +  33 

o  o 

Thus,  the  parameters  describing  the  burst  are  W,  (the  bubble  energy) 
Z,  (the  hydrostatic  head);  and  p  (the  fluid  density).  Dimensionally, 
these  are  singly  energy,  energy  density,  and  mass  density,  or  rather: 


V  has  dimensions  ML 


Z  has  dimensions  M 


p  has  dimensions  M 


(Energy) 

(Pressure) 

(Density) 


To  construct  a  "characteristic  length"  describing  the  explosion,  we 
may  combine: 

1 1 isf  Iff2]  (v/z)l/3  «  (the  characteristic  length,  feet) 

H1  /  “  I  (2'3) 

/  ml2  \  l/3(  yj}P16  I  m  ll/g  .  w1/3  £/s  __  t  (the  ctoracterlrtlc  tlM> 


‘Kl5/6 

H_V 

1  M  j 

Uh 

seconds) 


Other  characteristic  units  may  then  be  formulated,  for  example: 
Characteristic  area  =  Qp=(w/z)2^ 

Characteristic  volume  =  c?  =  w/z 
Characteristic  mass  =  po?  =  p  W/Z 
Characteristic  velocity  =  q/t  =  (z/p)1^2 


(2A) 


Characteristic  acceleration  =  a/ t2  =  (Z^^/W^^p) 
Characteristic  momentum  =  p  g^/t  =  W(p/Z)^2 
Characteristic  energy  =  p  Q?/t2  =  W 
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Characteristic  force  =  p  a* /V2  =  W2^ 

and  so  on. 

It  has  been  observed  that,  if  the  bubble  is  approximately  spherical 
(that  is,  gravity  and  interface  effects  are  small),  the  maximum  radius 
for  the  first  cycle  is  given  by 

R  =  12.6  a  (feet)  (2.5) 

and  the  first  oscillation  period  is: 

T  =  4.36  x  (seconds)  (2.6) 

It  might  now  be  pointed  out  that  W  may  be  more  simply  defined  as 
"that  weight  of  TNI  which  would  produce  the  'equivalent  bubble'  as 
regards  R  and  0/'  (or  a  and  r)* 

So  far,  we  have  considered  all  explosion  bubbles,  whether  point- 
source  or  high-explosive,  to  be  equivalent  after  the  time  of  shockwave 
separation.  This  is  a  fairly  good  assumption,  at  least  for  gross  be¬ 
havior  (the  non-homogeneity  of  the  steam  bubble  as  regards  temperature, 
density,  and  explosion  product  distribution  should,  nevertheless,  be 
kept  in  mind).  However,  the  fact  that  the  high-explosive  bubble  atmos¬ 
phere  is  composed  of  largely  non-condensable  reaction  products,  whereas 
the  point  source  bubble  contains  steam,  begins  to  have  its  effect  at 
about  the  time  of  the  first  bubble  minimum. 

This  "first  minimum"  is  accompanied  by  instabilities  of  the  bubble 
interface,  and  turbulence,  both  inside  and  outside  the  bubble.  In 
general,  energy  losses  occur  at  '  .Vble  minima  as  a  consequence  of  these 
phenomena.  In  the  high -explosive  bubble  case,  the  rise  in  the  internal 
pressure  causes  the  emission  of  a  pressure  pulse,  which  carries  energy 
away  from  the  system.  Also  (and  more  important)  the  turbulent,  unstable 
flow  near*  the  bubble  minimum  causes  transfer  of  heat  from  the  bubble 
atmosphere  to  the  water  environment  -  also,  the  kinetic  energy  involved 


in  this  turbulence  is  derived  from  the  bubble  energy  and  hence  repre¬ 
sents  an  additional  loss.  Consequently,  the  bubble  energy  available  in 
the  second  cycle  is  noticeably  less  than  that  of  the  first,  and  in  the 
HE  case,  is  about  39$  of  the  first  cycle  bubble  energy.  In  the  case 
of  a  steam  bubble,  all  the  above  losses  occur,  and  in  addition,  energy 
is  lost  through  another  mechanism.  As  the  bubble  collapses  and  the 
interface  becomes  unstable,  small  jets  of  water  penetrate  the  bubble, 
cooling  the  bubble  atmosphere  and  causing  condensation  of  steam,  thus 
depositing  the  latent  energy  of  condensation  in  the  water  as  heat.  This 
loss  is  quite  significant,  since  in  the  point-source  case,  the  bubble 
energy  available  for  the  second  cycle  is  only  about  8$  of  that  original¬ 
ly  available.1 

The  maximum  radius  and  the  period  af  oscillation  for  the  second 
and  subsequent  cycles  may  be  found  simply  by  replacing  the  total  bubble 
energy  (w)  by  the  reduced  amount  available  for  the  next  cycle.  Both 
cases  (HE  and  point-source)  are  compared  in  Fig.  2-1. 

THE  EFFECTS  OF  GRAVITY  AND  NEIGHBORING  SURFACES  -  BUBBLE  MIGRATION 

In  most  cases  of  practical  interest,  the  above  model  is  not  suf¬ 
ficient  to  completely  describe  the  phenomena.  So  fax,  we  have  assumed 
that  all  fixed  and  free-surfaces  are  far  away,  and  that  the  effect  of 
gravity  is  negligible.  An  extreme  case  of  the  failure  of  the  first  of 
these  assumptions  is,  for  example,  the  nuclear  test  Crossroads  Baker, 
in  which  the  shot  depth  was  a  great  deal  less  than  the  theoretical  first 
maximum  bubble  radius.  Also,  as  large  yields  become  of  interest,  even 
for  explosions  at  great  depths,  the  effects  of  gravity  are  certainly 
not  negligible. 

In  order  to  .  eal  with  the  effects  of  gravity  and  nearby  surfaces, 
various  models  have  been  advanced.  These  tend  to  divide  themselves 
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into  two  types,  which  will  be  dealt  with  separately. 

Spherical  Models 

Classically,  the  effects  of  gravity  and  surfaces  have  been  treated 
as  separable.  The  basic  assumption  involved  is  that  the  bubble  remains 
spherical,  and  that  gravity  and  surface  effects  may  be  treated  as  first- 
order  perturbations  of  the  non-migrating  case  discussed  above.  The 
effect  of  gravity  is,  of  course,  the  buoyancy  of  the  bubble.  Qualita¬ 
tively,  it  may  be  seen  that  the  total  amount  of  upward  momentum  generat¬ 
ed  by  the  bubble  through  its  first  cycle  due  to  gravity  is: 

T 

5  “  i*  *  buoyant 
o 

The  buoyant  force  is  given  by: 

F  =  4  *pg  (R(t))3 

3 

So  that,  T 

P  -  4  *pg  I  (R(t))3  dt 

where, 

p  =  upward  momentum 
T  is  the  first  bubble  period 
R(t)  is  the  bubble  radius  (a  function  of  time) 
p  is  the  fluid  density 

Now  since 

R  ~  OC 
T  ~  T 

we  have 
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If  ve  normalize  this  momentum  by  dividing  by  the  "characteristic 
momentum'  given  by  (po?)(o/ t),  where  a  and  t  are  cur  "characteristic 
lengtn"  anu  "characteristic  time",  we  see  that: 

^reduced  ~ 

That  is,  the  effect  of  gravity  increases  with  yield,  and  decreases  with 
shot  depth. 

Now,  this  momentum  is  accumulated  with  time  as  tbe  cycle  proceeds. 
When  the  bubble  is  large,  the  "virtual  mass"  cf  the  bubble  is  large, 
and  consequently,  the  rate  of  migration  is  small.  Near  the  bubble 
minima,  however,  the  virtual  mass  of  the  bubble  is  small,  and  conse¬ 
quently  the  upward  velocity  of  the  bubble  is  large. 

Superimposed  upon  this  buoyant  upward  migration  i3  the  effect  of 
adjacent  surfaces.  Intuitively,  this  is  somewhat  more  difficult  to 
grasp.  An  excellent  discussion  is  to  be  found  in  Ref.  8.  The  effect 
of  a  free  surface,  such  as  the  surface  of  the  water,  is  to  "repel"  the 
bubble  -  that  is,  to  cause  it  to  migrate  downward,  whereas  the  effect 
of  a  rigid  surface,  such  as  the  sea  bottom,  is  to  attract  the  bubble. 

Tbe  effects  of  eurfaces  increase  rapidly  as  the  surface  is  approach¬ 
ed.  Thus,  generally  speaking,  for  a  relatively  shallow  burst,  the  ef¬ 
fect  of  the  free  surface  will  predominate,  and  the  bubble  will  migrate 
downward.  As  shot  depth  is  increased,  the  "repulsion"  of  the  free  sur¬ 
face  drops  off  rapidly,  so  that  the  effect  of  gravity  begins  to  predomi¬ 
nate.  At  some  intermediate  shot  depth  called  the  "upper  rest  point" 
the  effects  cf  the  free  surface  and  the  buoyancy  of  the  bubble  cancel; 
the  bubble  does  not  migrate,  and  behaves  like  the  non-migrating  bubble 
described  previously.  Below  +his  depth,  migration  will  be  upward, 
until,  as  the  sea  bottom  is  approached,  the  "attraction"  of  the  sea 
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bottom  will  again  cancel  the  buoyancy:  this  is  the  "lower  rest  point". 
If  shot  depth  is  further  increased,  migration  will  be  downward. 

The  amount  of  migration  affects  the  amount  of  energy  lost  at  the 
bubble  minimum,  arid  consequently  it  has  been  necessary  to  evaluate  the 
loss  as  a  function  of  the  strength  of  migration.  This  has  been  done 
empirically  in  Ref.  7>  for  high-explosive  bubbles,  and  in  Ref.  1  for 
steam  bubbles. 


The  qualitative  picture  for  the  spherical  bubble  model  is,  then, 
as  follows.  The  bubble  is  assumed  to  remain  stationary  until  the  time 
of  the  first  bubble  maximum,  and  then,  as  the  bubble  contracts,  it  mi¬ 
grates  more  and  more  rapidly,  with  peak  velocity  being  attained  at  the 
time  of  the  bubble  minimum.  Upon  re -expansion,  the  bubble  slows  down, 
becoming  almost  stationary  at  the  time  of  the  second  maximum.  Thus, 
the  bubble  continues  to  migrate  in  a  series  of  "jerks",  remaining  spher¬ 
ical  throughout,  until  it  either  erupts  from  the  surface  or  dissipates 
all  its  energy. 


Quantitatively,  the  migration  between  the  first  and  second  bubble 

n 

maxima  may  be  represented  by  the  following  (considering  the  free  sur¬ 
face  only): 


where 

=  Shot  depth  (ft) 

Dg  =  Jepth  of  second  bubble  maximum 
=  Shot  depth  plus  33  feet 
R^  =  First  maximum  bubble  radius 
S  =  The  "surface  correction  coefficient" 


(2.7) 
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The  value  of  the  "surface  correction  coefficient"  is  a  natter  of 
sane  debate.  Estimates  vary  between  0.1  and  0.5.  Perhaps,  as  suggested 
in  Refs.  8,  9>  the  proper  value  of  S  is  a  function  of  the  shot  geometry. 
More  will  be  said  of  this  later. 


In  order  to  better  define  shot  geometry,  we  will  define  the 
"geometrically  scaled  depth",  as  follows: 


where  (Z  =  D  +  33) 


Note  that  (3  is  just  equal  to  the  shot  depth,  measured  in  first  maximum 
bubble  radii,  times  a  factor  of  12.6.  That  is,  a  charge  fired  at  a 
depth  such  that  £  =  12.6,  would  produce  a  bubble  that,  at  the  first 
maximum,  would  presumably  be  tangent  to  the  pre-shot  surface:  a  shot 
at  p  =  25  •  2  would  have  a  bubble  with  first  maximum  radius  equal  to  one- 
half  the  shot  depth,  and  so  on. 


Now,  equation  (2.7)  may  be  transformed  as  follows: 


Mg  =  specific  migration  -  ^1  -  ^2  =  12.4 


^\1/zh 
in]  ' 


S  |X2*6 j  Z 

!  e  I  B 


The  specific  migration  is  essentially  the  vertical  displacement  of 
the  bubble  center  between  the  times  of  the  first  and  second  bubble  max¬ 
ima  normalized  to  the  first  maximum  bubble  radius.  The  first  factor  in 
parenthesis  represents  the  gravity  migration,  and  the  second  is  the  so- 
called  ''free  surface  effect".  If  the  specific  migration  is  positive, 

the  bubble  migrates  upward  -  if  negative,  downward.  Fig,  2.2  presents 

-4-2  4 

specific  migration  as  a  function  of  0  for  W  =  10  ,10  ,1,  100,  10  , 

6  3 

10  and  10  pounds.  Note  that,  as  yield  increases,  the  effect  of 
gravity  gets  more  and  more  important.  It  should  be  pointed  out  that 
for  (3  <  ~15,  a  "res*-  point"  should  not  be  expected  even  if  predicted 
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by  this  theory,  since  the  surface  is  too  close  for  this  model  to  apply. 
In  order  to  generate  these  curves,  the  "surface  correction  coefficient" 
vas  assumed  to  be  a  function  of  £  which  gets  large  near  the  surface, 
but  approaches  0.1  at  large  depths  as  shown  in  Fig.  2-3*  This  informa¬ 
tion  was  derived  from  rest  point  data  from  various  sources,  as  noted. 

This  kind  of  "spherical  migration"  model  has  several  advantages 
and  has  many  useful  applications.  First,  it  predicts  with  reasonable 
accuracy  the  locations  and  times  of  emission  of  "bubble  pulses"  -  the 
pressure  pulses  emitted  at  bubble  minima.  The  results  derived  from 
this  model  agree  well  with  observations  when  migration  is  relatively 
weak,  that  is,  for  small  yields  fired  at  fair  distances  from  surfaces, 
or  for  large  yields  fired  at  great  depths.  Also,  it  is  analytic  and 
readily  applied,  requiring  no  more  than  pencil  and  paper  to  obtain 
results. 

On  the  other  hand,  there  are  many  shortcomings  to  this  sort  of 
model.  For  example,  it  breaks  down  badly  when  surfaces  are  nearby,  and 
when  migration  is  strong.  Clearly,  in  these  kinds  of  situations,  the 
assumption  of  a  spherical  form  (particularly  for  the  second  and  later 
maxima)  is  rather  poor.  Finally  (and  more  important,  at  least  as  re¬ 
gards  the  subject  matter  of  this  report )  the  model  says  nothing  what¬ 
ever  about  the  actual  events  occurring  within  the  bubble  near  tne  mini¬ 
mum. 


Non-Spherical  Models 

Essentially,  the  problems  involved  with  the  spherical  bubble 
models  arise  from  their  basic  assumption  that  the  bubble  retains  its 
spherical  form  throughout  all  phases  of  its  motion.  In  order  to  evalu¬ 
ate  the  effect  of  gravity  and  surfaces  upon  the  bubble  form,  models 
have  been  constructed  on  the  basis  of  the  following  set  of  assumptions: 
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1.  Spherical  symmetry  is  not  assumed,  but  social  symmetry  about 
the  shot  axis  normal  to  the  free  surface  plane  is  maintained. 

2.  All  flow  is  radial  about  the  shot  point. 

3.  The  position  of  the  bubble  interface  is  a  function  of  both 
time  and  0  (the  latitude). 

6,10 

There  are  two  such  models  in  use  at  present.  The  first 
starts  at  zero  time  with  a  sphere  of  gas,  the  size  of  the  initial 
bubble,  with  a  specified  internal  energy.  This  system  is  divided 
into  a  number  of  segments  with  various  values  of  0.  The  equations 
of  motion  are  then  solved  numerically  for  each  segment  (periodically 
re-evaluating  the  internal  energy  by  averaging  over  all  the  segments), 
keeping  in  mind  the  variability  of  the  hydrostatic  pressure  with 
depth.  The  results  are  as  follows: 

1.  The  bubble  is  non-spherical,  in  general,  and  an  assumption 
of  spherical  form  gets  worse  as  the  first  cycle  proceeds. 

2.  The  maximum  radius  and  the  "period"  are  functions  of  0.  In 
particular,  it  is  clear  that,  at  0  =  90°,  the  maximum  radius  and 
"period"  at  that  latitude  will  be  the  same  as  that  predicted  by  a 
spherical  model  for  the  bubble  as  a  whole,  that  is, 

R  (x/2)  =  12.6a 
T  (x/2)  =  4.36t 

However,  at,  say  0  =  l80°  (the  bubble  bottom),  the  period  and  radius 
will  be  reduced,  due  to  the  increase  in  hydrostatic  pressure  with 
depth.  In  the  upper  hemisphere  (say,  0  =  0°  -  the  bubble  top)  the 
maximum  radius  attained  will  in  general  be  greater  than  that  at  any 
other  value  of  0,  At  shallow  shot  depths,  however,  the  period  of  the 
bubble  top  may  be  quite  short,  due  to  the  small  mass  of  material  be¬ 
tween  the  bubble  and  the  water  surface.  At  great  shot  depths,  cm  the 
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other  hand,  this  effect  becomes  unimportant,  until  for  sufficiently 
deep  shots,  the  period  at  0  =  0°  may  be  the  longest  of  all,  due,  once 
again,  to  the  variation  of  hydrostatic  pressure  with  depth. 

The  second  non-spherical  model^  is  essentially  a  simplification 
of  the  first;  it  considers  only  the  bubble  top  and  bottom  (0  =  0°,l80°) 
and  neglects  the  internal  energy  of  the  bubble.  The  output  consists 
of  the  maximum  radii  attained  by  the  bubble  top  and  bottom  and  the 
times  at  which  these  maxima  occur.  The  results  are  essentially  ident¬ 
ical  to  those  of  the  more  general  model,  discussed  above.  Typical 
output  for  W  =  10^  pounds  is  presented  in  Fig.  2 -4a  and  2 -4b. 

Neither  of  these  models  deals  with  events  occurring  after  the 
first  minimum  -  in  fact,  the  models  terminate  prior  to  the  end  of  the 
first  cycle.  In  the  more  general  model,  termination  occurs  when  some 
segment  of  the  bubble  has  returned  to  the  explosion  point.  By  this 
time,  however,  the  important  information  about  bubble  behavior  has 
been  obtained.  The  fact  is  that  explosion  bubbles,  in  general,  col¬ 
lapse  asymmetrically.  For  shallow  bursts,  the  "top  period"  is  shorter 
than  the  "bottom  period",  whereas  for  deep  bursts,  the  reverse  is 
true.  This  asymmetric  collapse  gives  rise  to  a  jet  of  water  which 
passes  through  the  bubble^,^,^,^,^^,  thus  causing  translation  of  the 
system  as  a  whole. 

A  "characteristic  number"  describing  this  phenomenon,  correspond¬ 
ing  to  the  specific  migration  in  the  spherical  model,  is  the  "period 
ratio",  r^,,  where 

rT  =  Tl/2  ^  =  °)/Tl/2  ^  =  l8°0) 

That  is,  if  r^  <  1.0,  the  bubble  will  collapse  from  the  top,  and  we 
have  downward  migration;  if  r^,  >  1.0  we  have  "bottom  collapse"  and 
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the  system  translates  upward.  The  magnitude  of  should  also  cor¬ 
relate  with  the  strength  of  the  migration. 

Actually,  the  two  types  of  models  (spherical  and  non-spherical) 
are  describing  the  sane  phenomenon.  That  is,  for  deep  bursts,  the 
driving  force  is  bubble  buoyancy,  3nd  this  is  the  mechanism  for  the 
bottom  collapse  -  the  difference  between  the  top  and  bottom  periods 
in  the  non-spherical  model  is  a  consequence  of  the  pressure  gradient 
in  the  water,  and  buoyancy  is  simply  the  hydrostatic  pressure  inte¬ 
grated  around  the  surface  of  the  bubble.  Note  the  similarity  between 
Fig.  2-5,  which  shows  r^  as  a  function  of  (B  for  several  values  of  W, 
and  fig.  2-2,  which  has  ordinate  Mg.  When  the  spherical  model  indi¬ 
cates  appreciable  upward  migration,  r^  is  quite  a  bit  greater  than  1, 
whereas  when  downward  migration  occurs,  r^  is  less  than  1. 

The  most  obvious  shortcoming  of  these  non-spherical  models  is, 
of  course,  that  they  do  not  deal  with  later  phases  of  the  event  at 
all.  In  order  to  predict  events  after  the  first  recompression,  it  is 
necessary  to  resort  to  a  spherical  model. 

The  radial -flow  constraint,  while  an  improvement  over  a  spher¬ 
ical  constraint,  causes  problems  when  free  surfaces  are  nearby,  and 
also,  in  general,  just  prior  to  the  termination  of  the  model.  Under 
these  circumstances,  the  flow  has  a  tangential  component,  which  is 
not  accounted  for  in  the  model. ^  The  simplified  version  of  the 
model,  by  not  dealing  witu  times  after  maximum  expansion  and  by  treat¬ 
ing  only  the  bubble  top  and  bottom,  avoids  some  of  these  problems, 
but  the  imposition  of  a  radial -flow  velocity  potential  may  account 
for  the  model's  failure  to  predict  rest  points  and  other  observable 
data  precisely.  For  example,  unpublished  Hydra  series  film  data 
indicates  that  a  1-pound  charge  has  a  rest  point  at  8  feet,  but  the 
model  yields  r^  =  0.9  for  this  yield  and  depth,  which  would  indicate 
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downward  migration.  For  greater  scaled  depths,  however,  these  "free- 
surface"  problems  should  not  be  so  significant. 

Furthermore,  the  non-spherical  models,  in  contrast  to  the  spher¬ 
ical  models,  require  a  substantial  amount  of  computer  time  to  solve 
a  single  yield-depth  configuration.  If  these  models  are  further  re¬ 
fined  in  the  future,  however,  some  improvement  may  be  possible.  Never¬ 
theless,  the  elimination  of  the  spherical  constraint  represents  a  con¬ 
siderable  improvement  over  classical  models  from  a  conceptual  stand¬ 
point.  That  is,  the  radial -flow  bubble  model  provides  insight  into 
the  processes  actually  occurring  within  the  bubble,  rather  than 
simply  idealizing  the  system  as  a  pulsating,  translating  sphere. 
Although  not  perfect,  in  terms  of  explaining  observed  phenomena,  it 
is  a  much  better  model  when  surface  in  „eractxonB  or  strong  migration 
are  dominant  modes  of  behavior.  Intuitively,  it  is  much  more  satisfy¬ 
ing  also  -  it  explains  more  clearly  the  reason  for  the  so-called 
"surface  repulsion"  for  example,  and  suggests  reasons  for  the  observed 
non-spherical  form  of  the  second  and  later  maxima  when  the  migration 
is  strong. 

THE  SCALING  OF  EXPLOSION  BUBBLES  IN  THE  VACUUM  TEST  TANK 

Having  discussed  at  some  length  the  phenomenology  of  underwater 
bursts,  we  will  now  address  ourselves  to  the  problem  of  studying 
these  phenomena.  For  practical  reasons,  it  is  simply  not  feasible  to 
study  the  bubble  from  a  large  explosion.  As  we  have  seen,  however, 
small  explosions,  which  are  more  amenable  to  investigation,  produce 
bubbles  which  behave  differently  from  those  of  large  ones.  That  is, 
the  specific  migration  Mg,  for  example,  is  yield  dependent  and  never 
reaches  large  values  for  small  charges.  The  problem  is,  then,  to 
discover  what  environmental  parameters  must  be  varied  in  order  to 
reproduce  the  behavior  of  bubbles  from  large  explosions  in  the  labor- 
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story.  The  relevant  governing  equations  -are  derived  and  discussed  in 
appendix  A. 

Briefly,  it  may  he  shewn  that  the  behavior  of  the  bubble  (pulsa¬ 
tion,  migration,  etc)  from  a  large  explosion  may  be  reproduced  on  a 
small  scale  by  appropriate  adjustment  of  any  or  all  of  (the  air 
pressure  over  the  water,  in  atmospheres),  g(the  force  of  gravity),  or 
p  (the  density  of  the  fluid  medium).  The  effect  of  a  condensable  bubble 
atmosphere  does  not  scale  well  near  the  bubble  maxima,  but  error  may 
be  minimized  by  requiring  that  the  hydrostatic  pressure  at  the  shot 
point  be  sufficiently  great  so  that  the  bubble  interface  does  not 
"boil"  near  maxima.  Near  minima,  on  the  other  hand  (where  the  phase- 
change  process  becomes  important),  evaporative  and  condensive  effects 
should  scale  well.  The  effects  of  viscosity  and  compressibility  of 
the  fluid  are  not  properly  reproduced,  but  it  may  be  shown  that  they 
are  not  important  to  the  bubble  motion.  (See  Appendix  A). 


If  a  model  explosion  is  to  occur  underwater  in  a  laboratory 
environment,  in  a  fixed  test  tank,  the  only  environmental  parameters 
which  may  be  varied  are  the  yield  of  the  model  explosion  (Wm),  the 
shot  depth  (Dm),  the  water  temperature,  and  the  pressure  over  the  water 
(P^).  Hereafter,  any  difference  in  density  between,  say,  sea  water 
and  fresh  water  will  be  neglected.  Under  these  circumstances,  the 
scaling  equations  as  derived  in  appendix  A  become: 

w 

Yield  scale  factor  =  =  P,-^  (2*9) 

W  A 


Linear  scale  factor  = 


Time  scale  factor  = 


(2.10) 

(2.11) 
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vhere  is  in  atmospheres 

As  can  be  seen,  a  given  set  of  the  parameters  W^,  D^,  and 
define  a  unique  prototype,  of  yield  (W  /P^)  and  depth  (D^/P^). 

The  problem  is  that  there  is  a  lover  limit  on  PA;  the  vapor 
pressure  of  vater  even  at  the  freezing  point  is  0.00 6  atmospheres. 

In  general  practice,  even  this  limit  is  difficult  to  attain  -  a  more 
reasonable  value  might  be  0.010  atmospheres.  As  a  consequence,  a 

Q 

linear  scale  factor  of  100,  and  a  yield  scale  factor  of  10°,  repre¬ 
sent  the  upper  limit  on  the  scaling  capability  of  the  system.  This 
may  not  seem  a  serious  limitation,  and  indeed  it  is  not,  for  many 

applications.  In  a  laboratory  environment,  however,  it  is  not  practi- 

-4 

cal  to  use  model  yields  much  greater  than  about  10  lb.  Thus,  using 
vater  as  a  test  fluid,  the  largest  prototype  yield  which  may  be 
scaled  is  about  5  tons  of  TNT;  however,  if  simulation  of  the  behavior 
of  the  bubbles  formed  by  an  underwater  nuclear  detonation  is  of 
interest,  required  prototype  yields  fall  into  the  kiloton,  or  even 
the  low  megaton  ranges. 

In  order  to  simulate  the  migration  of  explosion  bubbles  caused 
by  detonations  of  this  size,  it  is  necessary  to  have  recourse  to  ap¬ 
proximations.  In  particular,  if  explosions  in  the  "very  deep"  cate¬ 
gory  are  of  interest  (P  >  ~  50),  the  effect  of  the  free  surface  is 
not  terribly  important.  Thus,  it  is  not  unreasonable  to  forfeit 
geometric  scaling  (equality  of  0)  in  favor  of  scaling  of  gravitational 
effects.  One  approach  is  as  follows.  If  the  model  and  prototype 
bursts  are  sufficiently  deep  that  the  effect  of  the  free  surface  is 
small,  then  scaling  of  bubble  behavior  will  be  accomplished  if  the 
Freuds  numbers  for  both  shots  are  equal  (Appendix  A,  equation  A.8). 
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This  is  the  so-called  "fourth -root  scaling  rule"  for  very  deep  bursts, 
that  is  often  applied  in  high-yield  high-explosive  field  tests1*^'12 
The  procedure  is  simply  to  compute  the  fourth-root  scaled  depth  for 
the  model  explosion,  and  then  to  assume  that  a  full-scale  explosion 
w:  \~u,  v!as  the  same  value  of  this  parameter  vill  generate  a  bubble  that 
•t  ■  in  the  same  way  as  the  model  bubble  does.  The  distance  and 
Ale  factors  are  then  simply  the  ratio  of  the  a’s  and  t's  for 
v  -odel  and  prototype. 


The  problem  that  may  arise  with  this  rather  simple  approach 
is  that  although  the  effect  of  the  free  surface  may  be  negligible  for 
the  prototype,  it  will  in  general  be  more  Important  for  the  model.  A 
better  approximation  to  actual  behavior  (particularly  if  free -surface 
effects  are  important  for  either  model  or  prototype)  is  to  determine 
the  prototype  shot  depth  for  which  bubble  behavior  will  be  the  same 
as  the  model  by  appeal  to  one  or  another  of  the  bubble  migration 
theories  discussed  previously. 


One  approach  might  be  to  specify  the  specific  migration.  Mg,  as 
a  scaling  criterion;  tnat  is,  scaling  is  accomplished  if: 
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The  same  procedure  may  be  followed  using  r^,  the  period  ratio,  as  a 
scaling  criterion.  These  approximations  are  not  precise,  for,  as  we 
have  seen,  neither  of  them  is  entirely  trustworthy  as  regards  the 
actual  effect  of  the  free  surface  on  the  bubble  motion.  Nevertheless, 
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they  should  not  lead  to  serious  discrepancies  for  deep  bursts.  The 
judgment  cf  the  user  must  be  applied  to  determine  vhich  method  is 
best  for  the  particular  configuration  of  interest. 


2k 


NRDL-2U-65 


Fig.  2-3  The  Surface  Correction  Coefficient,  S,  As  a 
Function  at'  Beta  (L/a)  • 
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Fig.  2 -4a  Maximum  Radii  for  the  Bubble  Top,  Bottom,  and  Side  far 
W  =  10,000  Pounds  as  Functions  of  Depth  of  Burst. 

Fig.  2 -4b  Bubble  Top,  Bottom,  and  Side  "Half  Periods"  for 
W  =  10,000  Pounds  as  .  auctions  cf  Depth  of  Burst. 
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Fig.  2-5  The  "Period  Ratio"  (r_)  as  a  Function  of  Beta  (D/a) 
For  W  =  1,  102,  104,  io6  and  108  Pounds. 


CHAPTER  3 


EXPERIMENTAL  APPROACH 


CONCEPT 

The  problem  involved  is  that  of  the  identification  and  evaluation 
of  the  processes  by  which  the  explosion  products  are  redistributed  in, 
end  ejected  from,  the  explosion  bubble  as  a  consequence  of  bubble  pul¬ 
sation  and  migration.  In  this  study,  point-source  explosions  are  of 
interest  since  the  results  will  eventually  be  applied  to  the  effects 
of  deep  underwater  nuclear  explosions.  Previous  work'*  has  indicated 
that,  at  the  first  bubble  maximum,  all  of  the  explosion  products  are 
contained  within  the  bubble,  and  that  most  of  them  reside  within  a  con¬ 
centrated  central  core.  The  objective  in  the  present  work,  then,  is 
to  examine  the  distribution  after  migration,  reccmpression  and  re¬ 
expansion  of  the  bubble  have  occurred.  It  is  apparent  that  the  extent 
and  nature  of  the  redistribution  of  explosion  products  will  depend  on 
the  behavior  of  the  bubble,  that  is,  the  strength  of  the  bubble  migra¬ 
tion.  Therefore,  the  method  adopted  was  as  follows. 

Explosion  tests  were  carried  out  in  the  NRDL  low-yield  underwater 
explosion  tank.  Three  scaled  shot  configurations  were  selected,  includ¬ 
ing  a  rest-point  condition,  a  condition  characterized  by  weak  upward 
bubble  migration,  and  a  strong  upward  migration  condition.  Scaling 
was  accomplished  by  reauction  of  the  air  pressure  within  the  tank, 
as  discussed  in  Chapter  2.  The  explosive  source  used  was  an  electri¬ 
cally  exploded  gold  wire  suspended  between  electrodes  in  the  center  of 
the  tank.  The  energy  densities  of  these  explosions  were  quite  high. 
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so  that  they  may  be  considered  point-source. 


The  investigation  was  carried  out  by  physically  sampling  the 
bubble  and  its  environment  at  the  time  of  the  second  bubble  maximum. 
The  samples  were  then  analyzed  for  gold  content  by  neutron  activation. 
In  order  to  avoid  perturbation  of  the  bubble  behavior  by  the  sampling 
system,  the  sampling  device  was  small  and  only  one  sample  was  taken 
for  each  explosion.  Consequently,  a  large  number  of  shots  were  fired 
at  each  condition  in  order  to  build  up  an  array  of  explosion  products 
concentration  data,  which  could  then  be  viewed  as  a  whole. 

THE  EXPERIMENTAL  FACILITY 

The  apparatus  used  in  this  study  was  basically  the  same  as  that 
used  in  the  work  reported  in  Ref.  5  and  much  of  the  basic  equipment 
has  been  described  at  some  length  in  Ref.  13 .  Therefore,  the  facility 
will  be  described  only  briefly  here,  with  emphasis  upon  those  features 
which  have  been  significantly  changed. 

The  Test  Tank 


The  explosion  tank  itself  is  essentially  an  upright  circular 
cylinder,  of  height  12  feet  and  inside  diameter  of  5  feet  (See  Fig. 
3-1).  Three  viewing  windows  are  provided,  each  60  inches  nigh  and 
18  inches  wide.  These  are  of  1-inch  tempered  glass,  and  are  mounted 
in  heavy  rectangular  steel  frames-  Tne  firing  electrodes  are  mounted 
in  ports  5  1-  feet  above  the  tank  bottom.  Other  ports  provide  access 
for  electrical  power,  pneumtaic  lines,  etc.  for  auxiliary  instrumenta¬ 
tion,  and  also  for  a  thermometer  for  measuring  water  temperature. 

Since  the  explosion  point  is  fixed  in  the  tank,  the  shot  depth  is 
varied  by  adjustment  of  the  water  level. 
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Access  to  the  tank  interior  is  through  the  front  viewing  window. 
Therefore,  it  is  necessary  to  lower  the  vater  level  in  the  tank  after 
each  shot  in  order  to  reset  the  sampler  and  replace  the  exploding 
wire.  This  window,  as  well  as  all  other  apertures,  is  sealed  to  main¬ 
tain  vacuum  when  the  pressure  is  lowered  prior  to  a  shot.  The  explo¬ 
sion  tank  walls  tend  to  collapse  inward  slightly  under  vacuum  -  the 
total  change  in  tank  diameter  being  about  l/l6"  under  full  vacuum. 

This  deflection  was  found  to  vary  linearly  with  the  pressure  differ¬ 
ential,  and  to  be  reproducible  within  10$. 

The  capacitors  used  to  store  the  explosion  energy  are  the  same  as 

TO 

were  used  previously.  J  These  are  six  General  Electric  Type  l4F?56,  50 
kilovolt  clam  shell  type,  low-inductance,  high  capacitance  discharge 
condensers  arranged  co-axially  in  parallel  just  outside  the  tank,  pro- 
viding  a  total  capacitance  of  5«76  pf.  J  Switching  is  accomplished 
by  a  three -ball  gap  which  is  pneumatically  operated  and  actuated  by  a 
solenoid  valve.  The  capacitor  bank  is  charged  by  a  Sorenson  Model 
2060-50  higf  oltage  DC  power  supply. 

Ordinary  tap  water  was  used  for  this  work,  as  its  background 
level  of  gold  is  far  below  the  range  of  interest.  An  auxiliary  hold¬ 
ing  tank,  located  outside  the  building,  was  used  to  store  the  water 
when  the  water  level  in  the  explosion  tank  was  to  be  lowered.  When 
not  is  use,  the  water  was  constantly  recirculated  through  a  sand 
filter  and  a  thermostatically  controlled  refrigeration  unit  in  order 
to  prevent  the  growth  of  algae  in  the  system  and  to  keep  the  vapor 
pressure  low. 

A  fresh  tank  of  water  was  introduced  periodically  (every  10-20 
shots),  in  order  to  avoid  buildup  of  the  gold  background  level.  The 
background  was  monitored  by  taking  samples  of  the  water  in  the  tank 
at  various  times  during  the  series. 
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The  Sampling  Device 


The  sampler  used  to  investigate  the  explosion  product  concentra¬ 
tion  was  similar'  to  the  device  used  in  earlier  work.'’  Essentially, 
it  consists  of  two  basic  parts  -  a  sampler  body  made  of  two  co-axial 
metal  tubes,  open  on  both  ends,  and  a  rod  with  three  attached  pistons; 
a  large  piston  on  one  end,  and  two  smaller  ODes  at  the  other  (Pigs. 
3-2,  3-3).  When  high  pressure  air  is  released  by  the  solenoid  valve 
(A)  it  passes  into  the  rear  chamber  (B)  and  exerts  force  on  the  large 
rear  piston  (c),  thus  causing  the  entire  piston  assembly  to  move  back¬ 
ward.  This'  draws  the  forward  pistpns  (D,E)  into  the  forward  part  of 
the  smaller  tube,  thus  trapping  the  material  inside.  The  motion  con¬ 
tinues  until  piston  (e)  strikes  the  bumper  ring  (F).  The  excess  pres¬ 
sure  in  the  rear  chamber  is  relieved  by  a  pin-hole  size  port  in  the 
side  of  the  chamber,  and  also  by  blow -by  past  the  rear  piston  (c). 

The  pistons  are  fitted  with  O-rings,  so  that  the  sample,  once 
taken,  is  vacuum  sealed  in  the  sampling  chamber.  In  practice,  leakage 
did  occasionally  occur  (about  20$  of  the  time)  but  it  is  quite  likely 
that  this  leakage  occurred  fairly  slowly.  In  general,  if  leakage  oc¬ 
curred,  the  flow  would  be  expected  to  be  into,  rather  than  out  of, 
the  sample  chamber;  unless  the  leakage  occurred  within  the  first 
second  or  so,  the  explosion  products  in  the  adjacent  water  would  be 
quite  dilute  because  of  mixing. 

The  volume  of  the  sample  chamber  is  1.05  cubic  centimeters,  this 
relatively  small  volume  being  necessary  to  avoid  excessive  perturba¬ 
tion  of  bubble  behavior.  This  is  less  than  0.1$  of  the  volume  of  the 
bubbles  under  study.  The  possibility  existed  that  the  "effective" 
sample  volume  may  vary  with  the  hydrostatic  pressure  because  of  the 
high  speed  of  the  sampler  pistons.  In  preliminary  tests,  the  sampler 
was  therefore  operated  underwater  at  hydrostatic  pressures  as  low  as 
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0.03  atmospheres  and  the  captured  water  weighed.  In  no  case  was  the 

■a 

derived  effective  sample  volume  significantly  different  from  1.05  cm'  , 

Monitoring  of  the  sample  operation  time  was  accomplished  in  two 
ways.  A  small  microswitch  (See  Fig.  3-2,  3-3)  is  mounted  in  the 
sampler  body  in  such  a  way  that  its  actuating  pin  rides  on  a  cam  that 
is  attached  to  the  piston  shaft.  Thus,  prior  to  operation,  the  switch 
is  open,  during  operation  it  is  closed,  and  when  the  sampling  opera¬ 
tion  is  completed,  the  switch  is  again  open.  A  6v  battery  is  con¬ 
nected  in  series  with  the  switch,  and  the  voltage  is  monitored  on  a 
Textronix  5^  5  A  oscilloscope.  Thus,  both  the  sample  time  and  sample 
duration  may  be  determined.  The  sample  duration  -  full  open  tc  full 
close  -  with  a  driving  pressure  of  100  psi,  was  about  4  msec*  which  is 
quite  short  compared  to  the  bubble  periods  of  interest.  The  sample 
time  -  th°  time  between  the  explosion  and  the  operation  of  the  sampler 
-  was  taken  as  the  time  midway  between  the  beginning  and  end  of  sampler 
operation.  San^le  times  determined  in  this  way  could  be  read  within 
about  1$.  The  sample  time  could  also  be  determined  from  the  high¬ 
speed  motion-picture  records.  When  the  sampler  closes,  the  excess 
pressure  in  the  rear  chamber  causes  a  large  bubble  of  air  to  be  ener¬ 
getically  expelled  from  the  relief  port  in  the  side  of  the  sampler 
body.  This  occurs  at  about  the  end  of  sampler  operation,  and  is 
clearly  visible  in  the  motion  pictures.  Sample  times  determined  in 
this  way  agree  within  2  milliseconds  with  those  measured  by  the  micro¬ 
switch  circuit. 

The  electrical  pulse  that  opens  the  solenoid  valve  and  initiates 
sampler  operation  was  obtained  from  the  sawtooth  output  of  a  Textronix 
535A  oscilloscope.  The  scope  was  triggered  by  the  explosion  discharge, 
and  the  sweep  delay  was  adjusted  so  as  to  obtain  the  desired  sample 
time.  Sample  times  could  be  reproduced  within  about  3  ms  in  this  way. 
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The  sampler  was  pre -positioned  in  the  tank  prior  to  the  shot  in 
such  a  way  that  the  sampler  axis  was  parallel  to  the  electrodes  and 
perpendicular  to  a  radius  joining  the  vertical  explosion  axis  and  the 
mid -point  of  the  sampling  volume »  This  was  done  for  two  reasons. 

First,  it  is  assumed  that  the  explosion  product  distribution  is  axially 
symmetric,  and  hence  this  orientation  would  minimize  uncertainties  due 
to  finite  sampler  size,  as  the  sample  chamber  diameter  is  less  than 
its  length.  Also,  this  approach  eliminates  ambiguity  concerning 
"effective"  sample  position,  since  all  motion  of  the  sample  pistons 
is  then  tangential.  Second,  since  there  is  no  fluid  circulation  around 
the  explosion  axis,  this  orientation  tends  to  eliminate  "masking"  of 
the  sample  by  the  front  piston. 

Photography 

High  speed  motion  pictures  were  taken  of  each  event  to  determine 
the  position  of  the  sample  relative  to  the  bubble,  and  to  make  direct 
measurements  of  bubble  phenomena.  Three  cameras  were  available,  and 
at  least  two  were  used  for  each  shot.  Pertinent  camera  information 
is  included  in  Table  3-1  and  Fig.  3-^* 

The  primary  function  of  cameras  1  and  2  was  to  provide  three- 
dimensional  coverage  of  the  sample  position  relative  to  the  bubble. 

Each  of  these  also  provided  additional  data  as  well.  Ideally,  camera 
1  should  have  been  located  further  from  the  tank,  but  this  was  prevent¬ 
ed  by  space  limitations.  Therefore,  the  camera  2  records  were  used 
for  precise  distance  measurements  (bubble  radii  and  the  like)  to  avoid 
parallax.  As  can  be  seen  in  Fig.  3-^»  from  the  camera  2  position, 
the  bubble  was  backlighted.  Consequently,  details  of  the  internal 
bubble  structure  could  be  discerned,  whereas  from  the  camera  1  posi¬ 
tion,  the  bubble  was  front-  and  sidelighted,  and  therefore  the  bubble 
appeared  opaque.  On  the  other  hand,  the  camera  1  records  proved 
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valuable  in  that  regions  of  high  explosion  product  concentration  were 
visible.  The  explosion  products  consist  principally  of  small  particles 
of  black  gold  oxide,  and  these  are  distinguishable  from  the  white 
background  and  from  the  many  tiny  bubbles  that  appear  during  the  ex¬ 
plosion  which  also  appear  white.  From  the  backlighted  camera  2  posi¬ 
tion,  both  these  small  bubbles  and  the  explosion  products  appear 
black,  and  cannot  be  resolved. 

Camera  3  was  not  used  for  all  the  shots.  The  nature  of  the  camera 
mechanism  and  the  graininess  of  the  Tri-X  film  used  resulted  in  poor 
photographic  resolution,  and  therefore  the  records  were  not  as  useful 
for  quantitative  measurements.  The  high  speed  at  which  the  camera 
was  operated,  however,  allowed  much  better  time  resolution  than  was 
obtainable  from  cameras  1  and  2. 

Distance  scaling  for  all  photographic  records  was  accomplished 
by  me'ms  of  marks  placed  on  the  viewing  windows.  These  were  located 
in  such  a  way  that  their  separations  corresponded  to  6-inch  intervals 
at  the  explosion  axis  from  the  vantage  points  of  the  various  cameras. 

Air  Pressure 


The  required  reduction  of  air  pressure  in  the  tank  was  obtained 
by  means  of  a  vacuum  pump,  and  was  measured  with  a  Wallace  and 
Tiernan  model  1*4-64  mercury  manometer  which  incorporates  a  barometer  for 
determination  of  ambient  pressure.  This  equipment  made  it  possible 
to  obtain  the  required  pressure  within  0.0001  atmospheres.  Generally 
speaking,  it  was  found  that  in  the  time  between  the  final  setting  of 
the  pressure  and  the  explosion  (about  one  minute)  the  change  in  pres¬ 
sure  due  to  leakage  and  de-gassing  of  the  water  was  never  more  than 
about  0.0004  atm.  This  pressure  uncertainty  corresponds  to  about  4  mm 
of  water,  which  is  about  the  same  as  the  uncertainty  in  total  shot 
depth  measurement. 
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Explosion  Yield  Determination 


To  determine  the  total  energy  ,/ield  of  the  explosion,  the  same 
system  was  used  as  in  the  work  reported  in  Ref.  13*  Essentially,  the 
method  involves  measurt  ment  of  the  total  released  energy  through  know¬ 
ledge  of  the  discharge  circuit  resistance  with  the  explosion  gap  short¬ 
ed,  the  inductance  and  capacitance  of  the  circuit,  the  voltage  on  the 
condensers,  ana  the  .Logarithmic  decrement  of  the  underdamped  sinusoi¬ 
dal  discharge  current  waveform.  This  current  waveform  was  determined 
by  integration  through  an  RC  network  of  the  voltage  induced  in  a  coil 

placed  adjacent  to  the  discharge  circuit,  and  the  resultant  voltage 

13 

was  monitored  on  a  Tektronix  551  oscilloscope. 

The  yield  measurements  were  attempted  for  each  explosion,  but 
instrumentation  problems  were  frequently  encountered  so  that  traces 
'•-.re  obtained  for  only  2/3  of  the  shots.  The  results  for  those  shots 
for  which  data  was  obtained  showed  wide  scatter  in  the  derived  yields 
for  identical  shot  condi-cions .  This  scatter  resulted  from  at  least 
three  sources.  First,  the  traces  indicated  that  polarization  may  have 
occurred  in  the  capacitor  bank  for  some  of  the  shots.  Also,  apparently 
some  shift  occurred  in  the  vertical  amplitude  calibration  of  the  oscil¬ 
loscope,  which  was  not  detected  until  sampling  was  nearly  complete. 
Finally  (and  most  important)  the  method  in  its  present  form  is  intrinsi¬ 
cally  impr-  cise  because  the  derived  yield  is  extremely  sensitive  to 
small  errors  in  measurement  of  the  values  of  the  current  peaks  in  the 
discharge  waveform;  these  can  be  determined  from  the  oscilloscope 
camera  records  orly  with  limited  accuracy. 

The  parameter  of  real  interest  in  this  study,  however,  was  not 
the  total  yield  (W* )  but  the  TliT  equivalent  bubble  yield  (W),  as  was 
discussed  in  chapter  2.  The  bubble  yield  is  readily  determined  by 
measurement  of  the  first  maximum  side  bubble  radius,  according  to  equa- 
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tion  2.5.  These  measurements  were  made  from  the  camera  2  films  for 
all  shots  which  provided  explosion  product  concentration  data.  The 
derived  yields  showed  much  less  scatter  tiian  the  total  yields  obtained 
from  the  discharge  decay  rates.  As  a  matter  of  academic  interest,  the 
value  of  k  (w/W;  equation  2.1)  was  calculated  for  each  shot,  and  the 
mean  value  was 

k  =  0.9  +  0.2 

This  may  be  compared  with  the  value  of  0.62  obtained  in  Ref.  14  for  a 
different  wire  geometry. 

Pres sure -Time  Measurements 


A  Crystal  Research  l/8"  tourmaline  gauge  was  installed  in  the 
tank  approximately  one-foot  directly  below  the  explosion  point.  In 
order  to  obtain  an  overall  pressure-time  history  of  the  explosion,  the 
output  was  recorded  on  a  Tektronix  53 5A  oscilloscope,  so  that  the  bubble 
pulses  provided  immediate -readout  estimates  of  the  bubble  periods  and 
hence  the  explosion  yield.  More  important,  the  same  output  was  record¬ 
ed  on  a  much  shorter  time  base  on  a  Tektronix  5^ 5 A  with  the  sweep  delay 
adjusted  to  obtain  a  detailed  record  of  the  first  bubble  pulse. 

Several  such  records  were  obi  ..ined,  and  results  are  presented  in  the 
next  section. 

The  t  . rmaline  gauge  was  removed  from  the  tank  about  half-way 
through  the  test  series,  since  an  adequate  number  of  bubble -pulse  wave- 
foi.s  had  be»n  obtained,  and  also  since  it  was  found  that  the  photographic 
records  yielded  bubble  periods  with  greater  precision. 
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experimental  PROCEDURE,  (typical  run) 


The  nature  of  the  experiment  involved  reproduction  of  essentially 
the  same  explosion  many  times.  Consequently,  it  was  quite  important 
thet  the  firing  conditions  for  each  shot  be  made  as  identical  as  pos¬ 
sible. 


The  firing  electrodes  were  installed  in  the  tank,  and  the  gap  set 
at  O.50  inches.  The  wire  segment  was  then  mounted  in  the  holes  in 
the  electrodes  in  such  a  way  that  no  stress  was  placed  on  the  wire 
when  the  tank  deformed  under  vacuum.  The  wire  simply  slid  freely 
further  into  the  electrode  to  accommodate  the  deflection.  (Fig.  4-l). 

The  sampler  was  first  disassembled  and  cleaned  with  a  dilute  HC1 
solution,  and  then  was  lubricated  with  a  silicone  grease  and  re¬ 
assembled.  It  was  installed  in  the  tank  in  the  "open"  position,  the 
tank  was  sealed,  and  the  water  level  was  raised.  The  electrodes,  ex¬ 
ploding  wire,  and  sampler  installed  in  the  tank  are  shown  in  Fig.  3-5- 
If  the  shot  conditions  required  particularly  low  air  pressures,  the 
tank  had  to  be  filled  very  slowly,  since  any  great  surface  turbulence 
would  have  caused  air  to  go  into  solution  in  the  water.  When  vacuum 
was  subsequently  drawn,  this  air  had  then  to  be  pumped  out  of  solu¬ 
tion,  which  might  take  several  hours  to  accomplish. 

The  water  level  was  adjusted  to  approximately  the  correct  shot 
depth,  and  vacuum  was  drawn.  Once  the  desired  vacuum  was  attained, 
the  water  level  was  re-adjusted  (the  deflection  in  the  tank  walls 
caused  the  water  to  rise  slightly).  The  system  was  allowed  to  stand 
for  about  15  minutes  to  allow  outgassing.  The  vacuum  was  .re-adjusted, 
and  the  shot  was  fired. 


40 


Immediately  after  the  shot,  the  vacuum  was  released,  the  water 
was  pumped  into  the  auxiliary  tank,  and  the  sampler  was  removed.  The 
outside  of  the  sampler  was  cleaned  with  HC1,  and  the  sample  was  re¬ 
leased  into  a  7  ml  polyethylene  vial  which  had  been  previously  weighed. 
The  sample  was  then  weighed,  and  the  sample  chamber  was  subsequently 
flushed  into  the  sample  vial  with  IN  HC1.  The  sample  vial  was  then 
heat-sealed  with  a  soldering  pencil  to  prevent  leakage. 

After  about  30  or  ^0  samples  had  been  collected,  they  were  ship¬ 
ped  to  General  Atomic,  San  Diego,  for  neutron  activation  analysis. 
Y/ithout  further  manipulation,  the  sample  vials  and  a  gold  comparator 

standard  were  placed  in  the  General  Atomic  TRIGA  Mark  I  nuclear  reactor 

12 

and  irradiate  '  for  30  minutes  at  a  thermal  neutron  flux  of  1.8  x  10 

n/cm  -sec.  ‘rue  samples  were  then  permitted  to  decay  for  ot  least  one 

day  to  eli  .inate  shortlived  interferences,  and  then  were  counted  on  a 

198 

multi-channel  gamma -ray  spectrometer  for  Au  ,  with  gamma  energy  = 

4ll  Kev  and  half-life  =2.7  days. 

Typically,  results  from  the  analysis  procedure  were  available 
about  a  month  after  sample  shipment,  and  consequently  the  sampling 
schedule  was  divided  into  only  two  separate  series.  The  first  was 
exploratory  in  nature,  being  designed  primarily  to  determine  the  dis¬ 
tribution  in  a  gross  fashion.  The  second  was  planned  after  the  results 
of  the  first  were  obtained,  to  provide  sufficient  data  to  make  quanti¬ 
tative  estimates  of  the  explosion  product  distributions  in  the  three 
shot  conditions  of  interest. 
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Fig.  3-^  The  Sqilosion  Tank  Facility  Showing 
Photographic  Layout. 


Fig.  3-5  General  View  Gf  Bqplosion  Tank 
Interior  Prior  to  Shot 
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CHAPTER  4 


DETAILS  OF  EXPERIMENTAL  CONDITIONS,  DATA  REDUCTION, 
AND  ANALYSIS 


THE  SHOT  CONDITIONS 

The  three  shot  configurations  under  study  were  selected  to  pro¬ 
vide  three  distinctly  different  sorts  of  bubble  behavior.  The  defin¬ 
ing  parameters  for  these  shot  conditions  are  summarized  in  Table  4-1. 

Since  the  independent  variable  involved  is  the  strength  of  bubble 
migration,  it  is  clear  that  a  shot  configuration  in  which  the  bubble 
does  not  migrate  at  all  is  of  interest.  As  was  discussed  in  Chapter  2, 
there  are,  in  practice,  two  ways  to  generate  a  non-migrating  explosion 
bubble.  First,  the  non-migrating  situation  is  approached  as  a  limit 
as  shot,  depth  is  increased  to  very  large  values  (assuming  that  the 
sea  bottom  is  always  very  far  away).  The  problem  involved  in  this  ap¬ 
proach,  however,  is  that  there  will  always  be  some  migration,  even 
though  it  may  be  extremely  weak,  and  that  very  great  shot  depths  are 
required.  This  means,  of  necessity,  that  the  bubble  will  be  quite 
small,  and  that  the  bubble  period  will  be  very  short.  In  the  experi¬ 
mental  situation,  sufficiently  great  shot  depths  are  simply  *.ot  avail¬ 
able  in  the  explosion  tank,  and,  in  any  case,  the  bubbles  would  be  far 
too  small  and  would  pulsate  too  rapidly  to  make  study  feasible. 

The  alternative  to  this  method  of  avoid5.ng  bubble  migration  is  to 
take  advantage  of  boundary  effects  -  in  particular,  the  "repulsion"  of 
the  free  surface.  An  explosion  fired  at  such  a  depth  that  the  surface 
correction  term  in  equation  2-7  is  zero,  i.e.. 
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will  generate  a  bubble  at  the  "upper  rest  point",  and  the  bubble  vill 
not  migrate.  Consequently,  the  following  shot  parameters  were  chosen 
for  Condition  A. 

The  air  pressure  over  the  water  was  set  at  0.1467  atmosphere,  and 
the  water  depth  was  adjusted  to  14.10  inches.  The  explosion  yield 
was  originally  intended  to  be  210  mg  TNT  and  was  determined  precisely 
a  posteriori  by  taking  the  mean  yield  of  all  data  shots  -  the  result 
was  that  the  TNT  equivalent  bubble  yield  was  209  +  12  mg  TNT.  These 
particular  shot  conditions  were  chosen  since  they  scale  the  upper  rest- 
point  of  a  one-pound  charge,  which  has  been  previously  documented. 

The  actual  prototype  conditions  obtained  were  W  -  0.994  +  0.057  lb, 
and  D  =  8.01  ft. 

The  second  shot  conf iguration  (Condition  B)  used  in  this  study 

was  chosen  to  be  characterized  by  upward  migration,  but  this  migration 

was  to  be  relatively  weak  (M  ~  0.8).  The  air  pressure  was  again  0.1467 

6 

atmospheres,  and  the  shot  depth  was  25. lo  inches.  The  explosion  yield 

was  agaii  intended  to  be  210  mg,  but  the  yield  obtained  by  averaging 

all  data  thots  of  this  series  was  173  +  11  mg  TOT  equivalent  bubble, 

so  that  the  specific  migration  (M  )  was  O.79.  This  disc,  ^pancy  in 

s 

yield  will  be  discussed  later  ?  i  oie  chapter.  It  was  hopci  that  this 

choice  of  shot  conditions  would  allow  correlation  of  results  with  the 

exploratory  work  done  at  the  second  bubble  maximum  reported  in  Ref.  5 

for  which  M  was  about  O.75.  More  will  be  said  about  this  later, 
s 

The  scaled  prototype  for  these  shot  conditions  is  a  charge  cf  bubble 
yield  0.325  +  0.05  pounds  fired  at  a  depth  cf  14.3  feet. 
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For  the  third  shot  configuration  (Condition  C),  it  vas  desired  to 
maximize  the  migration  strength,  and  consequently  the  air  pressure  *3S 
set  at  the  lowest  value  which  could  be  readily  reproduced,  to  obtain 
the  largest  possible  prototype  yield.  On  the  other  hard,  the  burst 
KPASt  be  kept  fairly  deep  (which  tends  to  repress  migration)  for  three 
reasons.  First,  as  shot  depth  is  decreased,  the  bubble  gets  quite 
large,  and  it  is  desirable  that  the  bubble  radius  remain  small  relative 
to  the  radius  of  the  explosion  tank  to  minimize  wall  effects.  Second, 
as  discussed  in  Chapter  2,  excessive  boiling  of  water  into  the  bubble 
may  occur  near  maxima  if  the  hydrostatic  head  at  the  shot  point  is  too 
lew.  Finally,  if  the  depth  of  burst  is  too  shallow,  since  the  migra¬ 
tion  is  upward,  the  second  bubble  maximum  may  be  perturbed  by  the  prox¬ 
imity  of  the  free  surface. 

The  shot  conditions  selected  for  condition  C  were  P^  =  0.0150 
atmospheres  and  D  =  29*00  inches.  The  designed  yield  was  138  mg  TNT. 
and  the  mean  bubble  yield  was  determined  to  be  134  +  13  mg  TNT.  The 
derived  specific  migration  was  1.47,  and  r^,  the  ratio  of  the  top  and 
bottom  bubble  half -periods,  was  1.20.  These  shot  conditions  correspond 
to  a  prototype  yield  of  5830  +  560  pounds  fired  at  a  depth  of  l6l  feet. 

Explosion  bubble  yields  for  individual  shots  were  determined,  as 
discussed  in  Chapter  3>  by  measurement  of  the  first  maximum  side  bubble 
radius.  The  TNT  equivalent  bubble  yield  is  then  given  by  Rg  =  12.6  a, 
op  in  this  case, 

W  =  0.01095  (R3)  (D  +  396  PA) 

where 

V  =  TNT  equivalent  bubble  yield  in  milligrams  TNT 
R  =  First  maximum  side  bubble  radius  (inches) 

D  =  Shot  depth  (inches) 

?A  =  Air  pressure  (atmospheres) 
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The  variation  from  shot  to  shot  in  the  total  hydrostatic  head  was 
less  than  0 which  is  quite  negligible  compared  to  the  variation  in 
yield.  Typically,  the  standard  deviation  in  yield  for  all  three  condi¬ 
tions  was  11-13  mg  TNT,  which  amounted  to  a  variation  of  about  &f> 
for  conditions  A  and  B,  and  about  10^  for  condition  C.  This  variation 
in  yield  wss  a  result  of  several  differenc  sources  of  e^ror.  First, 
the  firing  voltage  could  be  reproduced  only  within  about  +  500  volts. 
This  alone  would  tend  to  cause  a  variation  of  about  ^  6  mg  in  the 
bubble  yields.  Also,  the  three -ball-gap  trigger  resistance  depends  on 
the  breakdown  distance  in  air,  which  is  a  function  not  only  of  applied 
voltage,  but  also  of  ami  dent  air  density  and  absolute  humidity. 

Finally,  as  mentioned  in  Ref.  5,  the  length  of  the  exploding  wire  may 
have  varied  from  shot  to  shot,  but  it  is  felt  that  this  variation  was 
not  very  significant  in  the  present  work  -  one  standard  deviation  in 
wire  length  is  only  about  +  1 . 55b- 

Although  the  firing  voltages  for  conditions  A  and  B  were  identical 
(30  kv),  a  systematic  difference  between  the  yields  of  the  shots  was 
observed,  the  difference  between  the  mean  yields  being  about  three 
standard  deviations  of  the  individual  yields  for  either  shot  condition 
alone.  The  reason  for  this  difference  is  not  known,  since  the  only 
difference  in  the  nominal  shot  conditions  was  the  water  depth  in  the 
tank.  One  possibility  is  that  the  additional  mass  of  water  in  the 
tank  for  condition  B  may  have  caused  slight  mechanical  deformation  of 
the  tank  structure,  or  elastic  settling  of  the  tank  as  a  whole,  caus¬ 
ing  a  change  in  tor  geometry  of 'the  three -bail-gap  trigger,  or  even  a 
change  in  the  basic  circuit  parameters.  However,  an  understanding  of 
the  mechanism  for  the  variation  in  yield  is  not  critical  to  t  his  study, 
and  hence  the  matter  was  not  pursued  further. 
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DETERMUiATICH  OF  TOTAL  AMOUNT  OF  GOLD  PRESENT 


From  data  collected  concerning  the  effect  of  vacuum  cn  the  tank 
diametejj  it  is  known  that  the  gap  length  (separation  between  the  tips 
of  the  firing  electrodes)  was  0.446  +  0.006  inches  for  conditions  A 
and  B  and  0.438  +  0.006  inches  for  conditions  C.  The  exploding  gold 
wire  was  of  0.020  inch  diameter,  with  a  linear  density  of  99*6  mg/inch. 
Therefore,  if  the  total  mass  of  gold  involved  in  the  event  is  taken  as 
that  cf  e  segment  of  the  wire  of  length  equal  to  the  gap  length,  the 
"total  device"  or  total  input  gold  quantity  becomes  44.4  +  0.6  mg  for 
conditions  A  and  B  and  43.6  +  0.6  mg  for  conditions  C.  However,  there 
are  at  least  two  mechanisms  which  tend  to  change  this  amount.  First, 
the  "effective "  wire  length  is  somewhat  greater  then  the  gap  length  - 
erosion  of  the  electrode  tips  is  observed,  which  suggests  that  some 
wire  material  which  was  originally  within  the  electrode  was  involved 
in  the  explosion.  On  the  other  hand,  the  exterior  portion  of  the 
electrode  is  observed  to  be  coated  with  a  thin  film  of  material  after 
the  shot,  of  which  some  portion  is  undoubtedly  gold.  (See  Fig.  4-1). 

In  order  to  evaluate  the  net  effect  of  these  mechanisms,  a  seg¬ 
ment  of  wire  was  carefully  weighed  and  exploded  at  shot  conditions 
corresponding  to  condition  B.  The  electrode  tips  (containing  the  un¬ 
affected  ends  of  the  original  wire)  were  then  dissolved  in  aqua  regia 
and  analyzed  for  gold  content.  Tne  net  amount  of  gold  involved  in  the 
explosion  (the  total  initial  amount  less  that  present  in  the  tips)  was 
found  to  be  the  45.4  mg,  which  agrees  remarkably  well  with  the  above 
value.  This  indicates  that,  for  the  wire  geometry  under  study,  the 
effects  of  erosion  and  plating  tend  to  counteract  each  ether.  Hence, 
the  amount  of  gold  involved  in  the  explosion  for  all  three  conditions  was 
taken  as  45  mg,  and  this  value  is  probably  good  within  a  few  percent. 
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SAMPLING  DATA 


Sampling  Position 

The  position  of  the  sample  with  respect  to  the  bubble  was  determin¬ 
ed  from  the  camera  1  and  2  film  records*  Because  the  explosion  bubble 
motion  is  turbulent  and  the  initial  conditions  cannot  be  reproduced 
exactly,  the  bubble  shape  at  the  sec  end  maximum  varies  slightly  from 
shot  to  shot.  Consequently,  an  outline  of  the  ’mean*  second  bubble 
maximum  was  drawn,  and  then  the  position  of  each  sample  collected, 
with  respect  to  this  'mean1  bubble,  was  determined. 

The  vertical  coordinate  (2)  of  the  sample  with  respect  to  the 
actual  bubble  was  readily  found,  and  then  the  horizontal  cross-section 
of  the  actual  bubble  at  that  altitude  (z)  above  the  explosion  poir+- 
was  defined  by  the  four  interface  positions  provided  by  cameras  1  and 
2  -  this  cross-section  was,  in  general,  slightly  elliptical.  The  dis¬ 
tance  cf  the  sanqple  from  the  explosion  axis,  and  the  sample's  angular 
position,  were  then  determined  and  the  sample  placed  with  respect  to 
the  horizontal  bubble  cross-section.  Finally,  the  ratio  of  the  sample 
radius  (relative  to  the  center  of  the  cross-section)  to  the  bubble 
radius  at  that  altitude  and  angle  was  determined.  The  sample  position 
in  the  'mean*  bubble  was  defined  to  be  the  relevant  altitude  (z)  and 
a  radius  that  was  the  same  fraction  uf  the  'rean'  bubble  radius  at  z  as 
was  found  in  the  actual  bubble. 

Density  Measurements 

As  was  mentioned  in  Chapter  3 ?  the  mass  of  water  trapped  in  the 
soapier  was  measured  after  each  shot.  This  was  done  for  two  reasons. 
First,  the  amount  of  water  in  the  sample  determines  the  contribution 
of  the  gold  background  level  in  the  tank  to  the  gross  quantity  of 
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gold  in  the  sample  -  this  will  be  discussed  later.  Furthermore,  this 
measurement  is  valuable  in  itself,  as  it  provides  information  concern¬ 
ing  the  internal  structure  of  the  bubble. 

The  density  of  the  sample  vas  taken  simply  as  the  mass  of  the 
sample  divided  by  the  sample  chamber  volume  (l.05  cej).  Since  the 
bubble  atmosphere  is  5  *eam,  condensed  steam  would  tend  to  add  mass  to 
the  sample.  In  gene  1,  however,  the  mass  of  the  steam  is  negligible 
-  that  is,  if  the  internal  bubble  pressure  is  taken  as  0.2  atmos¬ 
pheres,  the  hydrostatic  pressure  at  condition  B  shot  depth  (which 
over-estimates  the  pressure  by  about  ar.  order  of  magnitude  or  so),  then 

the  mass  of  saturated  steam  in  the  sampler  volume  would  be  only  about 
_k 

10  grams,  which  is  about  the  lower  limit  of.  the  measurement  capa¬ 
bility  of  the  balance  used.  Consequently,  any  measurable  density 
detected  must  be  ascribed  to  the  presence  of  liquid  water,  or  spray. 

It  should  be  pointed  out  that  density  results  will  be  preferen¬ 
tially  low  since,  in  general,  it  is  not  possible  to  remove  all  of  the 
material  from  the  sample  chamber  without  flushing  the  sampler.  Also, 
since  spray  is  composed  of  droplets  of  finite  size,  the  results  will 
not  be  reproducible  -  a  droplet  of  radius  1  mm  weighs  4  mg,  and  a 
region  of  mean  density  0.1  grams/cm^  would  have  only  25  droplets  of 
this  size  per  ml.  Nevertheless,  in  terms  of  identifying  gross  phe¬ 
nomena,  these  measurements  proved  quite  valuable,  as  will  be  seen  later. 

Measurements  of  Explosion  Product  Concentration 

The  amount  of  gold  (in  micrograms)  present  in  each  sample  as 
reported  by  General  Atomic  was  initially  divided  by  1.05  cubic  centi¬ 
meters  (the  sample  chamber  volume)  to  obtain  gross  gold  concentra¬ 
tion.  Several  corrections  to  these  results  were  considered,  as  out¬ 
lined  below.  Relevant  data  for  all  samples  are  included  in  Tables 
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4-2,  4-3,  4-4,  and  4-5. 

The  Effect  of  Sample  Timing:  The  times  of  the  second  bubble 
maxima  were  39  msec,  37  msec,  and  220  msec  for  conditions  A,  B  and  C 
respectively,  and  the  time  of  sampler  operation  was  usually  vithln 
about  2  msec  of  these  maxima.  Of  course,  the  times  of  the  maxima  also 
varied  slightly,  due  to  the  variation  in  yield  from  shot  to  shot.  In 
general,  if  sample  time  was  appreciably  different  from  the  time  of 
the  second  maximum  for  a  particular  shot,  the  sample  was  not  used  for 
analysis.  Therefore,  it  is  felt  that  no  sample  timing  correction  to 
the  explosion  product  concentration  is  necessary. 

The  Effects  of  Yield  Variation:  The  appropriate  correction  factor 
to  be  applied  to  take  into  account  the  effect  of  shot-to-shot  varia¬ 
tion  in  yield  depends  on  the  mechanism  which  causes  the  variation. 

For  instance,  if  it  is  assumed  that  the  only  sources  of  yield  varia¬ 
bility  are  such  things  as  variation  in  firing  voltage  and  trigger 
resistance  -  that  is,  that  the  gap  length  is  constant  from  shot  to 
shot,  we  find  a  "characteristic  concentration" 

C  =  M/o? 

where  M  is  the  mass  of  the  exploding  wire  and  c?  is  w/z.  We  assume 
that  Z  and  M  are  constant,  so 

C'  =  c(w/w* ) 

where  C  snd  W  are  measured  values  and  C'snd  W’  are,  respectively, 
the  corrected  gold  concentration  and  the  mean  bubble  yield  for  the 
shot  condition.  If  we  assume,  on  the  other  hand,  that  all  variation 
in  yield  is  due  to  gap  length  variation,  then,  approximately 
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B-10  5-70  169  Pull  -  1.39  1-32  3-7  2.3 
B-ll  5 ‘71  170  57*8  0.06  3.89  3.70  2.7  2.8 
B-12  5.98  195  739-3  O.70  34.4  32.8  l.o  1.5 
B-13  5.90  187  53^.6  0.51  26.7  25. ^  1.0  2.0 


TABLE  4-4 

Sampling  Bata  -  Condition  C 


Shot 

Firat  Max 

Side  Bub. 
Bad.  (in.) 

V 

(ng  TMT) 

Saaple 

Water 

Content 

(=e) 

Sample 

Density 

(«/«*) 

Uncorrected 

Gold 

Qmntit> 

(us) 

Corrected 

Gold 

Concentration 

(ns/c=3) 

Location 
z(ln. )  r* in. j 

C-l 

6.94 

128 

575-7 

0.55 

0.632 

0.600 

17.6 

0.0 

C-2 

6.58 

109 

607.2 

O.58 

1.67 

1.59 

16.7 

0.0 

c-3 

6.85 

123 

JUll 

- 

1.73 

1.64 

16.0 

2.0 

C-4 

7.01 

132 

301.5 

0.29 

1.18 

1.12 

15-5 

0.0 

C-5 

6.89 

125 

489.1 

0.47 

2.64 

2.50 

15.5 

2.0 

c  -6 

7.03 

133 

551.8 

0.53 

0.0005 

0.000 

15-5 

4.5 

C-7 

7.07 

135 

247-9 

0.24 

0.806 

0.767 

15.0 

3.4 

C-8 

7-12 

137 

249.2 

0.24 

0.859 

0.817 

14.5 

0.0 

c-9 

7-00. 

132 

341.9 

0.33 

1.28 

1-22 

14.5 

3.2 

C-10 

6.98 

130 

156.6* 

0.15* 

1  37 

1.30 

13-7 

0.0 

c-u 

6.65 

112 

269-5 

0.26 

1-33 

1.27 

13.0 

0.0 

C-12 

7.01 

732 

184.9 

0.18 

0.681 

0.648 

13.O 

2.1 

C-13 

7.18 

141 

407.1 

O.39 

0.188 

0.178 

13.0 

5-1 

C-l4 

7  -12 

137 

-15-9 

O.U 

0.496 

0.471 

12.3 

3-5 

C-15 

7-07 

135 

8.9 

0.26 

0.516 

0.490 

12-5 

5.0 

c-l6 

7.00 

131 

-05-4 

0.27 

0.754 

0-717 

11.0 

0.8 

C-17 

7.0 

136 

522.0 

0.50 

1.17 

1.10 

11.0 

2.0 

C-l8 

7-2<* 

145 

67-2 

0.06 

0.312 

0.296 

11.0 

3.0 

C-19 

7-ce 

132 

71.3 

0.07 

0.365 

0.347 

11.0 

5.0 

C-20 

7.26 

146 

606.1 

0.53 

0.534 

0.507 

11.0 

5.8 

C-21 

7.01 

132 

461.7 

0.44* 

0.0177 

0.0151 

11.0 

7-0 

C-22 

6.6i* 

112 

Full 

- 

1.69 

1.6l 

9.5 

0.0 

C-23 

6.64 

112 

246.8 

0.23 

1.59 

1.51 

9.5 

0.0 

C-2** 

6.70 

115 

425-7 

0.4l 

0.120 

9.5 

1.4 

C-25 

6.97 

130 

103.0 

0.10 

0o75 

0.357 

9.5 

2.9 

C-2o 

7.12 

137 

Full 

- 

3.33 

3.21 

8.4 

0.0 

C -2/ 

7.20 

143 

Full 

- 

0.874 

0.827 

8.4 

4.0 

C-28 

6.76 

118 

Full 

- 

0.629 

0.595 

8.2 

5-7 

C-29 

6.97 

130 

310.6 

0.30 

1.89 

1.80 

7.0 

0.0 

c-30 

7.23 

147 

55-5 

0.05 

O.05O 

0.905 

7.0 

4.2 

c-31 

■7.76 

146 

599-2 

0.57* 

1-38 

1.31 

4-5 

0.0 

C-32 

6.95 

128 

718.0 

0.68* 

1.11 

1.06 

4.5 

2.0 

c-33 

7.32 

150 

1032.9 

0.98 

0.498 

0.466 

2.0 

0.0 

c-31* 

7-33 

151 

Full 

1.0 

0.465 

0.438 

2.0 

3.0 

« 

Lovor  Unit 
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TABLE  4-5 

Late -Time  Condition  C  Sampling  Data 
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and  we  find  that  no  correction  is  necessary. 

The  fact  is  that  variation  in  gap  length  is  responsible  for  some 
fraction  of  the  yield  variation,  but  the  importance  of  gap  length  vari¬ 
ation  relative  to  other  sources  has  not  been  quantitatively  determin¬ 
ed.  In  any  case,  the  variation  in  yield  itself*  was  small  (at  worst, 
about  +  lOjo  for  condition  C),  and  consequently  the  correction,  if 
applied,  would  be  even  smaller.  Therefore,  no  correction  of  concen¬ 
trations  for  individual  shot  yield  was  used. 

The  Effect  of  Background:  The  choice  of  gold  as  the  exploding 
wire  material  was  made  for  two  reasons.  First,  quantities  of  gold 
may  be  determined  by  neutron  activation  analysis  with  great  precision 
even  in  trace  concentrations.  Second,  as  a  consequence  of  gold's  rel¬ 
ative  scarcity  in  nature,  the  problems  of  contamination  in  sample 
handling  and  background  are  minimized. 

In  order  to  prevent  the  buildup  of  a  high  background  level  in 
the  tank  due  to  deposition  in  the  water  of  the  explosion  products 
after  each  shot,  a  new  tank  of  water  was  periodically  introduced,  as 
discussed  in  Chapter  3*  If  it  is  assumed  that  all  the  explosion  pro¬ 
duces  are  uniformly  mixed  after  each  shot,  in  the  total  amount  of 
water  in  the  system  (1500  gallons ) ,  then  the  upper  limit  on  the  rate 
of  background  buildup  is  about  8  nanograms  of  gold  per  cubic  centi¬ 
meter  per  shot.  Consequently,  in  the  first  sampling  series,  a  7  ml 
sample  of  the  water  was  taken  after  introduction  of  the  fresh  tank 
before  any  shots  were  fired,  and  another  such  sample  was  taken  just 
prior  to  emptying  the  tank.  It  was  found  that  the  background  level 
of  gold  in  fresh  tap  water  was  quite  low  and  fairly  constant  -  O.OO37 
pg/ml  +  31$.  Therefore,  for  the  second  shot  series,  no  samples  were 


60 


taken  of  the  fresh  tanks;  the  concentration  vas  assumed  to  be 
initially  about  0.004  pg/ml. 

The  intended  procedure  was  to  assume  a  linear  buildup  of  the 
background  for  each  tank  of  water  of  the  form 


Cbg'°B0  +fln 
o 

where  C2  is  the  cone .ntration  in  fresh  tank  water,  and  n  is  the 
I3Go 

number  of  shots  fired  in  that  tank  of  water.  The  value  of  A  (the 
buildup  rate)  is  obtained  from  the  sample  taken  just  prior  to  empty¬ 
ing  the  contaminated  tank  of  water.  Thus,  the  background  concentra¬ 
tion  of  gold  in  the  water  is  defined  for  each  shot,  and  consequently 
the  contribution  of  this  background  to  the  gross  amount  of  gold  in 
the  sample  is  just  the  background  concentration  times  the  volume  of 
water  found  ia  the  sample  chamber. 


Upon  examination  of  the  background  data,  however,  it  became  ap¬ 
parent  that  the  background  concentration  was  being  affected  not  only 
by  the  rate  at  which  explosions  were  contributing  gold  to  the  water, 
but  also  by  mechanisms  that  were  extracting  the  material.  The  gold 
explosion  products  were  not  in  solution,  but  were  in  the  form  of  an 
extremely  finely  divided  particulate.  Continuous  recirculation  of 
the  water  through  the  filter  over  long  periods  of  time  apparently  ex¬ 
tracted  a  significant  amount  of  material  from  the  water;  gold  particles 
may  also  have  adhered  to  irregularities  in  the  pipe  walls,  and  so  on. 
Since  sufficient  data  were  not  available  to  permit  evaluation  of  the 
effects  of  these  transfer  processes,  the  background  concentrations 
derived  from  the  linear  fits  to  the  data  discussed  above  were  used, 
with  some  reservation. 


Fortunately,  the  contribution  of  the  background  to  the  total 
amount  of  gold  in  most  of  the  samples  was  negligible.  If  the  fraction 
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of  the  gross  concentration  represented  by  background  is  computed  for 
the  samples,  the  following  table  results: 

Percentage  of  Samples 

Fractional  Background  With  Fractional  Background  Contri- 

Contribution  (percent)  but ion  Less  than  that  Value. 

0.01  9.5 

0.1  38.1 

1  84.1 

10  95.2 

The  three  samples  with  the  background  concentrations  equal  to  or 
greater  than  10$  of  the  total  were  samples  C-21  (10$),  B-l  (l4$),  and 
C-6  (330$)*  Fortunately,  C-21  and  B-l  were  each  fired  as  the  first 
shot  in  a  fresh  tank  of  water,  and  the  background  is  well  known  under 
these  conditions.  As  can  be  seen,  in  Table  4-4,  shot  C-6  provided  a 
sample  containing  less  gold  than  would  be  predicted  for  a  sample  of 
water  containing  only  the  calculated  background  level.  The  only 
reasonable  conclusion  is  that  the  region  in  which  this  sample  was 
taken  was  virtually  uncontaminated,  except  for  background. 

Sources  of  Error  in  Explosion  Product  Concentration  Data:  Errors 
in  the  analysis  of  the  samples  were  of  great  concern,  since  Ref.  5 
reports  variations  of  as  much  as  a  factor  of  2  in  standard  samples 
sent  to  check  precision  of  the  neutron  activation  analysis.  In  order 
to  minimize  these  errors,  the  sumple  handling  procedure  was  simplified, 
as  described  in  Chapter  3,  and  it  appears  that  this  solved  most  of  the 
problems.  General  Atomics  Inc.  reported  standard  deviations  in  their 
reported  values  of  gold  content  which  were  generally  under  +  1$. 

These  deviations  were  determined  from  counting  statistics  only,  how- 
ver.  In  order  to  evaluate  overall  analysis  errors,  such  as  variations 
in  neutron  flux,  etc.,  a  series  of  standard  samples  were  prepared  and 
included  with  test  samples,  the  gold  contents  varying  over  the  range 
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0.01  pg  -  80  pg.  The  reported  values  were  within  10$  of  the  correct 
values  for  these  standards,  indicating  that  analysis  problems  need  be 
of  little  concern. 


Another  source  of  error  and  non-reproducibility  is  intrinsic  in 
the  explosion  products  themselves.  The  explosion  products  are  not  ell 
in  solution,  bu:  are,  at  least  in  part,  in  the  form  of  small  gold 
oxide  particles.  This  introduces  the  possibility  of  "particle  count¬ 
ing  error"  -  that  is,  at  low  levels  of  concentration,  the  sample  may 
contain  only  a  few  particles,  and  the  addition  or  subtraction  of  even 
one  particle  may  make  a  noticeable  difference  in  th<»  drived  gold 
concentration.  In  general,  with  a  uniform  particle  ...ize,  the  uncer¬ 
tainly  (l  a)  associated  with  a  measurement  of  n  particles  is  Vn.  If 
we  assume  that  all  the  explosion  products  are  associated  with  identi¬ 


cal  spherical  gold  oxide  particles  of  diameter  D  microns,  that  the 
density  of  the  particles  of  Au„0  is  3*6  grams/cm,  and  that  the 
sample  volume  is  1.05  cm,  then  the  density  of  the  gold  in  the  parti- 
cles  is  3.2  grams/cm,  and  it  may  be  shown  that  the  percentage  uncer¬ 


tainty  in  the  gold  concentration  is: 


~  X  lOOjt  =  +  0.13 


where  the  concentration  is  in  micrograms  per  cubic  centimeter. 


As  can  be  seen,  the  uncertainty  is  strongly  dependent  on  particle 
size.  To  investigate  the  problem,  a  sample  was  taken  of  the  particles, 
the  sample  was  passed  through  a  10  mp  millipore  filter,  and  the  par¬ 
ticles  were  photographed  (see  Fig.  4-2).  Virtually  all  of  the  particles 
are  below  1  p  in  size.  Even  if  1  p  is  tsken  as  the  mean  particle  size, 
the  uncertainty  over  the  concentration  range  of  interest  is  under  +  1$, 
and  consequently  the  problem  of  particle  size  is  not  relevant. 
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Another  source  of  error  arises  from  the  use  of  the  sampling 
technique.  First#  the  sample  position  must  he  determined,  and  it  is 
estimated  that  this  may  he  done  vithin  0.1  inch.  Second,  although  the 
sampler  size  is  small  relative  to  the  huhble,  the  sample  is  of  finite 
extent  -  the  diameter  of  the  sample  chamber  is  0.435  inch.  This  tends 
to  cause  an  ’’averaging”  of  the  gold  concentration  data  across  the  region 
in  vhich  the  sample  was  taken.  This  effect  will  not  he  important  ex¬ 
cept  where  concentration  gradients  are  very  large,  and  may  indeed 
have  the  beneficial  effect  of  masking  turbulent  microstructure  which 
would  otherwise  tend  to  obscure  the  general  distribution.  Neverthe¬ 
less,  it  may  be  safely  stated  that  the  average  concentration  found  in 
the  sample  will  reflect  the  exact  concentration  at  sane  point  along 
the  sandier  diameter,  in  all  likelihood  fairly  close  to  the  center. 
Somewhat  arbitrarily  therefore,  the  total  uncertainty  in  "effective" 
sample  position,  reflecting  both  uncertainty  in  actual  position  and 
the  effect  of  finite  size  was  taken  as  +  0.2  inches. 

It  should  be  pointed  out  that  intrinsic  uncertainty  is  intro¬ 
duced  by  the  nature  of  the  experiment.  That  is,  the  hydrodynamic  pro¬ 
cesses  occurring  in  the  underwater  explosion  bubble  phenomenon  are 
turbulent  in  nature,  and  consequently  small  differences  in  initial 
conditions  tend  to  generate  larger  differences  from  shot  to  shot  as 
the  event  proceeds.  The  effects  of  the  inability  to  reproduce  the 
phenomenon  exactly  are  extremely  difficult  to  estimate  a  priori. 

Ideally  under  these  circumstances  the  best  procedure  would  be  to  ac¬ 
cumulate  many  data  for  each  sailing  position  within  the  bubble,  and 
then  to  construct  an  a  posteriori  estimate  of  uncertainty.  Quite 
aside  from  the  difficulty  of  reproducing  a  sample  position,  however, 
this  procedure  would  involve  the  taking  of  such  a  large  number  of 
samples  as  to  render  the  project  impossible,  from  the  standpoints  of 
both  cost  and  time.  The  choice  involved  was  essentially  between  hav¬ 
ing  relatively  few  data  positions  for  which  the  uncertainty  is  well 
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known  quantitatively,  or  having  a  larger  number  of  data  with  the  un¬ 
certainty  not  so  well  defined;  this  investigator  chose  the  latter  ap¬ 
proach,  and  the  results  appear  to  justify  this  choice.  At  one  of  the 
sample  positions,  where  two  samples  (C-22  and  C-23)  were  taken,  th2 
gold  concentration  results  were  1.6l  and  1.51  (igrams/cnr  respectively. 
Although  this  close  agreement  may  be  somewhat  fortuitous,  it  does  in¬ 
dicate  that  the  uncertainty  caused  by  irreproducibiiity  is  probably 
no  more  important  than  the  total  caused  by  the  sources  of  error  al¬ 
ready  discussed.  Furthermore,  although  the  rigorous  procedure  of  ob¬ 
taining  numerous  data  for  each  sampling  position  was  not  followed,  a 
large  number  of  samples  was  obtained  at  various  positions,  such  that 
the  sample  density  in  the  bubble  was  fairly  high.  It  seems  reason¬ 
able  to  assume  that  the  effect  of  irreproducitility  is  not  terribly 
important  if  the  sample  data  give  rise  to  meaningful,  continuous  con¬ 
centration  gradients  and  contours,  which  indeed  they  do.  The  judgment 
of  the  user  must  then  he  employed  to  evaluate  the  significance  cf 
deviations  in  the  data.  It  should  also  be  pointed  out  that  too  large 
a  value  at  one  sample  position  may  be  accompanied  by  too  small  a  value 
at  an  adjacent  position,  so  that  when  the  concentration  data  is  inte¬ 
grated  over  many  sample  positions  to  obtain  a  total  gold  content  in 
some  region,  the  error  in  the  total  will  in  general  be  less  than  the 
error  in  the  individual  sample  data. 
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Fig.  4-2  Explosion  Product  Particles  From  The 
Submerged  Exploding  Wire 
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CHAPTER  5 

RESULTS  AND  CONCLUSIONS 

DISTRIBUTION  AND  TRANSPORT  MECHANISMS 

The  problem  involved  in  this  study  is  to  determine  how  the  central 
distribution  of  explosion  products  found  in  Ref.  5  for  the  first  bubble 
maximum  is  changed  by  the  hydrodynamic  mechanisms  of  bubble  collapse, 
re-expansion,  and  migration.  For  the  initial  discussion,  we  will  con¬ 
sider  the  three  shot  conditions  separately. 

The  Non-Migrating  Case  -  Condition  A 

The  Distribution  a"  Explosion  Products:  Ideally,  a  non-migrating 
explosion  bubble  should  remain  spherical  even  after  several  cycles. 

The  condition  A  bubbles.,  however,  at  the  second  maximum  took  the  form 
of  an  ellipsoid  of  revolution  centered  at  the  shot  point  with  the 
vertical  radius  (~  3*5  inches)  being  greater  than  the  horizontal 
radius  (,>2.5  inches)  as  shown  in  Fig.  5-1*  This. sort  of  behavior  has 
been  observed  also  at  the  prototype  shot  conditions  (Unpublished  HYDRA  I 
film  data),  and  is  probably  a  consequence  of  the  fact  that  the  non- 
migrating  situation  was  obtained  by  "balancing"  the  effects  of  the  free 
surface  and  gravity  rather  them  by  firing  at  so  great  a  shot  depth 
that  the  pressure  differential  across  the  vertical  bubble  diameter  was 
negligible  relative  to  the  total  head.  It  was  hoped  that  this  depar¬ 
ture  from  spherical  foim  would  not  seriously  affect  the  distribution. 

Sampling  was  attempted  along  both  the  vertical  and  the  horizontal 
bubble  radii.  Unfortunately,  it  was  found  that  sample  positions  above 


or  below  the  shot  point  would  generally  induce  upward  or  downward 
migration  toward  the  sampler.  After  numerous  attempts,  one  sample 
(A -3)  was  obtained  along  the  vortical  radius  of  a  non-migrating  bubble. 
If  sample  radii  were  normalized  to  the  bubble  radius  at  the  appropriate 
latitude,  it  was  found  that  this  sample  concentration  agreed  quite 
nicely  with  data  taken  along  the  horizontal  radius.  It  was  therefore 
concluded  that  the  contours  of  equal  concentration  were  ellipsoids  of 
revolution  of  the  same  eccentricity  as  the  bubble  itself.  The  sample 
positions  and  corrected  concentration  data  are  shown  in  Fig.  5-2,  and 
the  reduced  radial  concentration  is  presented  in  Fig.  5-3* 

In  order  to  more  easily  compare  results  of  shots  at  one  condition 
to  those  at  another,  as  well  as  to  simplify  integration  for  total  con¬ 
tent,  the  position  and  concentration  data  may  be  readily  transformed  to 
"reduced"  values.  That  is,  we  may  take  a  (i.e.  (W.Z)1^)  as  our  basic 
unit  of  length;  the  unit  of  volume  then  becomes  a? (i.e.  W/z).  Rather 
than  use  units  of  mass  (in  this  case,  jig  of  gold)  as  a  measure  of  the 
total  quantity,  we  may  use  simply  fractions  of  the  "total  device",  or, 
as  a  basic  unit  of  concentration  tid/oP  (micro-devices  per  cubic  alpha). 
If  these  transformations  are  made,  using  1  d  =  45  mg  gold  and  a  =  0.510 
inches,  we  obtain: 

1  ng/cm^  =  48.49  nd/cP 

Nov,  we  may  refer  our  ellipsoidal  hubhi°  to  a  spherical  bubble  of 
the  same  volume.  If  this  is  done,  we  obtain  that  the  equivalent 
spherical  radius  at  the  second  maximum  is: 

R2  =  5-48  a 

(Recall  that  the  radius  at  the  first  maximum  is  12.6a).  Earlier  work 
concerning  the  distribution  of  explosion  products  within  the  bubble  at 
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the  first  maximum  indicated  that  at  this  time,  most  of  the  explosion 
products  reside  in  a  concentrated  central  core."*  This  result  corrobo¬ 
rates  earlier  theoretical  work  which  indicated  the  same  sort  of  behavior 
for  deep  underwater  nuclear  explosion  bubble s.^  The  study  of  the  first 
tabble  maximum  was  done  in  a  manner  similar  to  that  of  the  present 
wo~k,  using  an  exploding  g'  Id  wire,  physically  sampling  the  bubble, 
and  analyzing  for  gold  by  neutron  activation.  The  central  concentra¬ 
tion  was  found  to  be  more  than  three  orders  of  magnitude  higher  than 
the  concentration  near  the  periphery,  the  concentration  away  from  the 
center  dropping  off  approximately  with  l/r^.  Therefore,  the  distri¬ 
butions  in  reduced  units  for  the  first  and  second  bubble  maxima  are 
as  in  Fig.  5-^«  If  this  distribution  is  integrated,  the  fraction  of 
the  total  device  lying  within  a  given  reduced  radius  as  a  function  of 
reduced  radius  may  be  obtained  -  this  may  then  be  compared  with  the 
situation  at  the  first  bubble  maximum,  as  shown  in  Fig.  5-5  • 

As  can  be  seen,  whereas  all  of  the  explosion  products  reside 
within  the  bubble  at  the  first  bubble  maximum  for  the  non-migrating 
case, only  about  18$  are  inside  the  bubble  at  the  second  maximum. 

Most  of  the  remainder  lie  within  a  fairly  short  distance  from  the 
interface,  but  are  outside  the  bubble.  Just  adjacent  to  the  inter¬ 
lace,  concentrations  are  sufficiently  high  so  that  the  explosion  pro¬ 
ducts  are  clearly  visible  on  the  camera  1  film  records.  The  con¬ 
taminated  "patches"  at  the  interface  are  shown  in  Fig.  5-6. 

The  Mechanisms  of  Explosion  Product  Transport:  As  was  discussed 
in  Chapter  2,  as  the  explosion  bubble  collapses  from  its  first  maxi¬ 
mum  expansion,  the  internal  pressure  rises  rapidly.  Shortly  before 
the  bubble  minimum,  the  gaseous  bubble  atmosphere  begins  to  change 
phase  near  the  interface  -  mass  is  lost  to  the  liquid  medium.  Essen¬ 
tially,  the  interface  begins  to  move  with  respect  to  coordinates 
fixed  in  the  fluid.  Consequently,  the  explosion  products  which  were 
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embedded  In  the  atmosphere  are  transferred  through  the  bubble  inter¬ 
face  into  the  liquid  environment.  This  transfer  is  visible  on  the 
c  amera  1  film  records  as  a  "blackening"  of  the  originally  white  bubble 
interface  just  prior  to  the  minimum. 

As  the  contraction  proceeds  to  the  minimum,  more  and  more  material 

is  lost.  It  should  be  pointed  out  that  the  minimum  bubble  radius  is 

rather  poorly  defined,  since  the  internal  pressure  and  temperature 

are  well  above  the  critical  point  for  water  at  and  near  the  minimum. 

Also,  at  the  bubble  minimum  the  bubble  is  irregular  and  "spiky"  in  ap- 

9 

pearance  due  to  Taylor  instability.  This  irregularity  is  reflected 
in  the  "blotchiness"  of  the  explosion  products  visible  near  the  inter¬ 
face  at  the  second  maximum. 

At  the  bubble  minimum,  the  quantity  of  explosion  products  within 
the  bubble  may  be  very  small.  As  the  bubble  re  '-“oands,  however, 
the  adjacent  water  begins  to  vaporize  due  to  the  u.jp  in  internal  pres¬ 
sure,  and  material  is  thus  transported  back  into  the  bubble;  the  inter¬ 
face  is  now  moving  outward  with  respect  to  the  fluid.  Most  of  the 
explosion  product  material  remains  outside  the  bubble,  however,  even 
after  the  phase  change  process  is  cccplete.  The  total  bubble  energy 
remaining  in  the  bubble  for  the  second  cycle  is  a  great  deal  less 
(~  8$)  than  that  which  was  available  in  the  first,  and  consequently 
the  total  amount  of  water  re-vaporized  upon  the  second  expansion  is 
appreciably  less  than  that  which  was  condensed  prior  t-o  the  minimum. 
Also,  near  the  bubble  minimum,  considerable  turbulence  is  generated, 
resulting  in  both  internal  and  external  mixing.  The  result  is  that, 
at  the  second  bubble  maximum,  there  still  remains  a  high  central  con¬ 
centration,  but  concentration  gradients  are  not  nearly  so  large  as  at 
the  first  maximum.  The  bulk  of  the  explosion  products  (~  82$)  now 
reside  just  outside  the  bubble  interface.  These  transport  mechanisms 
are  illustrated  in  Fig.  5-7  for  Condition  A. 
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Weak  Migration  -  Condition  B 


The  Explosion  Product  Distribution:  In  many  respects,  condition  B 
may  be  treated  as  a  perturbation  of  condition  A  (the  non-migrating 
case).  At  the  second  maxima  .  the  bubble  is  roughly  pear-shaped  as 
shown  in  Fig.  5-3*  The  pro*  .oe ranee  at  the  top  is  a  consequence  of 
the  bottom  collapse  which  occurs  just  prior  to  the  time  of  maximum  re- 
compression,  as  will  be  discussed  later. 

The  sampling  data  is  presented  in  Fig.  5-9*  In  order  to  convert 
this  data  to  "reduced'"  units  (as  was  done  for  condition  A)  the  appro¬ 
priate  conversion  factors  are: 

a  =  O.456  inches 
1  pg/cnf*  =  34.51  pd/o? 

The  concentration  along  the  vertical  centerline  axis  within  the  bubble 
is  shown  in  Fig.  5-10.  If  the  off-axis  concentrations  are  plotted  as 
a  function  of  (cylindrical)  sample  radius  normalized  to  cylindrical 
local  bubble  radius,  all  the  values  below  z/a  =  15  and  z/a  =  5  fall 
on  a  smooth  curve,  as  shown  in  Jig.  5-11.  A  few  relevant  values  ob¬ 
tained  from  Ref.  5  for  a  similar  second  maximum  (  see  Chapter  4)  also 
fall  on  the  same  curve.  The  similarity  in  form  of  Fig.  5-11  for 
condition  B  and  Fig.  5-3  for  condition  A  should  be  noted. 

The  region  above  z/a  =  15  (see  Fig.  5-9)  any  be  considered  as 
relatively  homogeneous  in  concentration  -  the  sampler  size  is  fairly 
large  relative  to  this  protuberance,  so  that  the  data  would  tend  to  be 
"average"  values.  If  the  mean  concentration  is  taken  as  200  pd/a^, 
the  total  content  within  the  bubble  below  z/a  =  15  is  4-6$  of  the 
device,  leaving  93  to  95$  unaccounted  for.  The  location  of  this  ex¬ 
cluded  fraction  is  not  difficult  to  determine,  however.  Although  no 
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sampling  data  is  available  outside  the  bubble  in  the  "blotchy"  region 
along  the  sides  visible  in  Fig.  5-12,  this  blotchiness  is  in  all  like¬ 
lihood  similar  to  that  observed  in  condition  A,  and  denotes  a  high  con¬ 
centration  of  explosion  products  just  outside  the  interface.  Even 
more  important,  the  heavily  contaminated  region  visible  in  Fig.  5-12 
just  below  the  bubble  was  sampled  at  two  positions}  the  mean  concen¬ 
tration  in  this  region  is  well  above  1000  pd/o?.  Although  sufficient 
data  is  not  available  to  obtain  a  detailed  distribution  in  these 
regions,  concentrations  of  this  magnitude  indicate  that  the  remainder 
of  the  device  is  to  be  found  there.  The  distribution  of  the  explo¬ 
sion  products  in  the  various  zones  of  interest  is  summarized  in 
Fig.  5-13. 

The  Mechanisms  of  Explosion  Froduct  Redistribution:  The  differ¬ 
ences  between  conditions  A  and  B  arise  from  the  collapse  of  the  bottom 
of  the  bubble  near  the  first  recompression  and  the  translation  of  the 
system  as  a  whole  in  the  latter  case.  As  has  been  seen  (Chapter  2), 
just  prior  to  the  recompression,  after  much  of  the  explosion  product 
material  has  been  lost  to  the  environment  due  to  phase  changes,  the 
bubble  bottom  passes  through  the  explosion  point  and  impinges  upon  the 
top.  The  upward  jet  of  water  thus  formed  scavenges  some  of  the  con¬ 
tained  explosion  products  and  penetrates  the  water  above  the  bubble. 

At  this  time,  the  gaseous  portion  of  the  bubble  is  translating  upwards 
with  high  velocity  through  the  water,  leaving  behind  it  much  of  the 
explosion  products  which  were  lost  due  to  phase  change  of  the  bubble 
atmosphere.  Thus,  as  the  bubble  re-expands  and  adjacent  water  is 
vaporized,  much  of  the  water  that  passes  into  the  bubble  atmosphere 
is  uncontaminated  and  hence  returns  no  explosion  products  to  the 
bubble.  Early  in  the  re-expansion  the  narrow  central  vertical  water 
jet  breaks  up  into  spray >  Is  dispersed,  and  is  in  all  likelihood  in 
large  part  vaporised,  since  at  the  second  maximum  no  central 
coluiur  was  detected  by  the  sample  density  measurements.  As  a  conse- 
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quence,  the  internal  distribution  at  the  second  maximum  is  similar  to 
that  in  the  non -migrating  ca3e  except  as  follows:  The  central  core 
has  been  smeared  vertically  through  the  bubble  by  the  upward  water 
jet,  moving  the  region  of  highest  internal  concentrations  up  into  the 
residual  protuberance  at  the  top.  The  total  internal  content  of  explo¬ 
sion  products  is  only  about  l/3  that  found  in  the  non-migrating  case, 
since  the  system,  upon  re-expanding,  has  translated  into  relatively 
uncontaminated  water.  Thus,  although  some  of  the  lost  material  is  to 
be  found  in  a  degenerate  form  of  the  external  irregular  "shell"  of 
condition  A  around  the  lower  half  of  the  bubble,  much  of  the  material 
remains  in  a  relatively  small,  highly  concentrated  patch  below  the 
bubble,  centered  at  about  the  position  at  which  the  bubble  minimum 
occurred.  These  transport  mechanisms  for  condition  B  are  illustrated 
in  Pig.  5-1^. 

Strong  Migration  -  Condition  C 

Much  more  data  was  accumulated  at  this  shot  condition  than  at 
any  other,  for  two  reasons.  First,  the  bubble  at  the  second  maximum 
is  physically  much  larger  than  was  true  for  condition  A  and  B,  and 
the  flow  is  more  complicated,  requiring  more  samples  for  adequate 
definition  of  the  distribution.  Second,  the  bubble  behavior  (migration 
strength)  is  characteristic  of  nuclear  explosions  fired  at  nominal 
depths  o?  burst. 

The  motion  picture  sequence  in  Fig.  5-15  shows  general  features 
of  the  bubble  phenomena;  relevant  sample  information  at  the  second 
maxinum  is  presented  in  Fig.  5-16.  The  data  obtained  concerning  the 
internal  density  is  shown  in  Fig.  5-17  and  indicates  the  existence  of 
a  centra^  column  or  jet  of  extremely  heavy  spray  (density  ~  0.25 
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grams /car)  which  passes  vertically  through  the  bubble.  This  column  cf 
spray  is  also  faintly  visible  on  the  backlighted  camera  2  film  records 


in  Fig.  5-15*  Clearly,  the  bubble  may  no  longer  be  considered  as  a 
gaseous  sphere,  but  is  actually  toroidal  in  form. 


As  was  done  previously,  position  ana  concentration  data  may  be 
transformed  to  reduce  quantities  by  application  of  the  conversion 
factors: 


a  =  O.56O  inches 
1  pg/cm3  =  63.76  pd/a? 

The  data  may  then  be  plotted,  and  contours  drawn  as  shown  in  Fig. 

5-18.  Several  features  of  the  explosion  product  distribution  may  now 
be  distinguished  as  illustrated  in  Fig.  5-19*  Near  the  top  of  the 
bubble  and  extending  slightly  into  the  adjacent  water  above  is  a 
region  (i)  of  relatively  high  explosion  product  concentration  (~  100 
pd/a3)  which  is  more  or  less  homogeneous.  This  region  corresponds  to 
the  top  of  the  central  spray  jet  shown  in  Fig.  5-17*  The  total  con¬ 
tent  of  this  region  is  about  7#  of  the  device.  Continuous  with  region 
I  is  a  vertical  hollow  cylinder  of  high  concentration  which  corresponds 
in  position  with  the  edge  of  the  central  water  jet,  extending  all  the 
way  to  the  bottom  of  the  bubble  (region  II ).  The  concentration  in 
this  region  varies  between  60  and  90  pd/a3,  and  the  total  contained 
quantity  of  explosion  product  is  about  3$  of  the  total  device.  Out¬ 
side  regions  I  and  II  but  within  the  bubble  is  a  torus-shaped  zone 
(region  III)  wherein  the  spray  content  is  quite  low,  and  where  the 
explosion  product  concentration  varies  from  about  30-50  pd/a3  around 
the  perimeter  to  about  20  pd/o3  at  the  center.  This  region  accounts 
for  about  10$  of  the  total  amount  of  explosion  products  present. 

Below  region  I  and  inside  region  II  is  a  relatively  thin  shell  (region 
IV)  where  the  concentration  is  lower  than  in  the  adjacent  water,  vary¬ 
ing  from  ~  50  pd/o?  near  the  top  to  values  below  10  pd/a3  near  the 
lower  bubble  interface.  The  total  amount  of  material  present  in  this 
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region  is  only  about  0.4$  of  the  device.  Contained  within  region  IV 
and  continuous  with  the  black  "skirt"  of  debris  visible  below  the 
bubble  in  Fig.  5-20,  is  region  V,  where  concentrations  are  again  high, 
reaching  a  peak  of  200  pd/a?  near  the  lower  bubble  interface,  and  de¬ 
clining  slowly  below  the  bubble.  Due  to  the  small  volume  involved  in 
that  portion  actually  within  the  visible  bubble  interface  (region  Va) 
the  total  content  involved  there  is  only  about  3$  of  the  total  device, 
but  in  the  region  below  the  bubble  (Vb)  the  total  volume  is  quite 
large  due  to  the  spreading  of  the  region  as  the  explosion  point  is 
approached.  It  is  therefore  assumed  that  the  remaining  amount  of 
the  device  (77$)  is  distributed  within  this  region.  Although  in¬ 
sufficient  data  is  available  to  integrate  a  total  content,  the  large 
volume  involved  and  the  relatively  high  concentrations  present  certain¬ 
ly  indicate  that  this  conclusion  is  not  unreasonable. 

The  Redistribution  Mechanisms:  In  condition  C  the  bottom  col¬ 
lapse  and  upward  translation  are  dominant  modes  of  behavior.  The  up¬ 
ward  jet  of  water  passes  through  the  bubble  and  strikes  the  top  at  a 
time  when  the  internal  bubble  pressure  is  not  yet  particularly  high. 
Consequently,  only  a  small  fraction  of  the  explosion  products  have 
been  lost  to  the  environment  when  the  jet  rushes  through  the  highly 
contaminated  center  of  the  bubble.  This  effect  can  be  seen  in  the 
pressure-time  records  shown  in  Fig.  5-21  taken  near  the  bubble  minima. 
The  "bubble  pulse"  from  condition  A  shots  is  merely  a  single  sharp 
peak,  whereas  for  conditions  B  euid  C  a  second,  lower  peak  is  observed 
just  prior  to  the  main  (gas  recompression)  peak.  This  early  pulse  is 
a  consequence  of  the  impact  of  the  upward  jet  against  the  upper  bubble 
interface,  and  has  also  been  observed  at  full-scale  field  tests'*''*. 

The  time -separation  between  the  peaks  was  about  0.15  msec  for  condi¬ 
tion  B  and  about  1.2  msec  for  condition  C,  reflecting  the  earlier  col¬ 
lapse  at  the  strong  migration  shot  configuration. 
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Therefore,  a  substantial  portion  of  the  explosion  products  becomes 
mixed  into  the  top  and  sides  of  the  jet  prior  to  the  time  of  maximum 
rec expression.  Physically,  the  "visible"  bubble  minimum  is  relatively 
large,  but  most  of  its  volume  consists  of  the  water  in  the  upward  jet, 
the  gaseous  portion  of  the  bubble  being  a  relatively  small  annulus. 
Events  occurring  near  the  minimum  are  shown  in  Fig.  5-22*  Once  again, 
as  the  bubble  recompresses,  most  of  the  remaining  material  is  lost  to 
the  environment,  and  is  swept  into  a  skirt-shaped  wake  behind  the 
bubble  as  it  translates  upward.  The  motion  after  the  first  recompres- 
sion  is  observed  to  be  somewhat  similar  to  that  of  a  ring,  or  spherical 
vortex,  with  the  strong  upward  flow  of  the  central  jet  continuing 
and  spreading  out  along  the  top  of  the  bubble.  This  circulating  flow 
causes  mixing  of  explosion  products  from  the  jet  into  the  toroidal 
gaseous  part  of  the  bubble,  and,  furthermore,  causes  a  small  portion 
of  the  material  which  was  initially  lost  near  the  minimum  to  be  drawn 
back  up  into  the  center.  In  contrast  to  the  behavior  in  condition  B, 
the  Jet  does  not  degenerate,  but  remains  a  dominant  feature  through¬ 
out  the  motion. 

DISCUSSION 

In  general,  as  should  be  expected,  the  distribution  at  the  second 
bubble  maximum  Is  more  dispersed  than  that  of  the  first.  Furthermore, 
as  migration  strength  is  increased,  the  explosion  products  are  not 
only  translated  upwards  to  a  greater  and  greater  extent  but  are  more  ani 
more  widely  distributed,  as  can  be  seen  in  Fig.  5-23,  where  the  dis¬ 
tribution  at  the  second  maximum  for  all  three  conditions  is  compared 
with  that  at  the  first  bubble  maximum.  The  total  amount  retained  in 
the  bubble  at  the  second  maximum  is  fairly  large  when  no  migration 
occurs,  but  drops  off  as  migration  strength  is  increased,  due  to 
bubble  re-expansion  into  uncontaminated  water.  As  migration  strength 
is  further  increased,  however,  the  upward  jet  becomes  a  dominant 


77 


feature  of  the  system,  and  the  retained  explosion  product  fraction 
again  increases.  Furthermore,  as  migration  strength  increases,  the 
kinetic  energy  of  migration  becomes  a  significant  term  in  the  bubble's 
energy  budget.  Thus,  the  internal  energy  at  the  minimum  decreases, 
recompression  becomes  less  extreme,  and  less  material  is  lost  through 
the  phase  change  process. 

Two  earlier  studies  have  been  made  concerning  the  distribution 

of  explosion  products  from  a  submerged  exploding  wire.  One  of  thes 

concerned  itself  with  the  distribution  at  the  first  bubble  maximum, 

and  has  already  been  discussed. ^  More  recently,  work  was  undertaken 

at  Malaker  Laboratories,  Inc.,  in  which  an  electrically  exploded  silver 

ribbon  submerged  in  light  mineral  oil  was  used  as  an  explosive 
l6  17 

source.  ’  1  The  system  could  be  evacuated  to  obtain  various  migra¬ 
tion  strengths.  Essentially,  the  procedure  involved  was  that  of  tak¬ 
ing  high-speed  motion  pictures  of  the  explosion  and  attempting  to 
track  the  "debris  cloud."  It  is  suggested  in  Ref.  17  that  the  debris 
resided  in  a  thin  shell  near  the  interface  at  the  first  bubble  maxi- 
rau  .  It  is  also  concluded  that  little  or  no  e::plosion  product  debris 
is  transported  by  the  bubble  and  released  only  after  oscillation  has 
ceased.  These  conclusions  are  at  variance,  both  with  the  present 
study  and  with  the  work  at  the  1st  bubble  maximum  presented  in  Ref.  5. 
There  are  several  possible  reasons  for  this  disagreement,  as  follows. 

First,  the  energies  involved  in  the  explosions  were  extremely  low- 
in  the  range  0.1  -  1  calorie  or  about  0.1  -  1  mg  THT  (as  compared  to 
^  150  mg  THT  for  the  present  work).  At  these  low  yields,  the  frac¬ 
tion  of  the  energy  required  merely  to  vaporize  the  wire  material  may 
well  become  a  substantial  fraction  of  the  total.  The  explosion  will, 
in  any  case,  be  of  fairly  low  energy  density  so  that  the  "point- 

source"  approximation  may  not  hold,  and  the  explosion  product  particles 

1  ^ 

may  be  fairly  large.  These  problems  would  in  all  likelihood  serious- 
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ly  affecb  the  initial  distribution  of  explosion  products  in  the  bubble. 
Also,  the  explosion  bubbles  were  strongly  perturbed  by  the  small  size 
of  the  container  relative  to  the  bubble  size,  and  the  large  firing 
electrodes,  as  shown  in  the  photographs  of  the  bubbles  in  Refer¬ 
ences  16  and  17. 

Furthermore,  the  use  of  photography  as  a  method  of  quantitative 
determination  of  explosion  product  distribution  can  be  misleading. 

This  is  particularly  true  if  backlighting  is  used  (see  Chapter  3). 

If,  for  example,  the  explosion  products  are  fairly  dilute  over  a 
large  region,  they  may  be  invisible  on  film,  but  the  integrated  total 
explosion  product  content  over  the  region  may  still  be  large.  On 
the  other  hand,  if  ejqplosion  products  are  concentrated  in  a  shell 
outside  the  interface,  they  may  be  indistinguishable  from  the  inter¬ 
face  itself. 

In  order  to  apply  the  present  results  to  possible  full-scale 
situations,  one  or  another  of  the  scaling  criteria  discussed  in 
Chapter  2  may  be  applied.  Strictly  speaking,  the  prototype  shot 
yields  and  depths  for  conditions  A,  B,  and  C  respectively  are  0-99  lb 
at  8.01  feet,  O.83  lb  at  14.3  feet  and  5830  lb  at  l6l  feet.  As  was 
pointed  out  in  Chapter  2,  however,  approximate  scaling  may  be  accom¬ 
plished  by  resorting  to  one  or  another  of  the  bubble  models. 

In  free-field,  condition  A,  the  non-migrating  case,  does  not 
correspond  to  any  full  scale  shot  conditions  for  large  yields;  it  is 
only  approached  as  a  limit  as  shot  depth  is  increased  to  very  large 
values,  except  for  yields  of  less  than  100  lb  or  so,  for  which  an 
upper  rest  point  is  possible.  If  the  sea  bottom  is  nearby,  a  lower 
rest  point  may  be  attainable  and,  if  the  bubble  is  not  excessively 
distorted  by  proximity  to  the  bottom,  condition  A  may  be  a  good  model 
for  bubble  behavior  and  explosion  product  distribution. 
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Conditions  B  and  C  may  be  referred  to  other  yields  and  shot 
depths  by  means  of  the  bubble  models.  That  is,  condition  B  is  char¬ 
acterized  by  M.  =  0.79  (that  is,  AD.  _  =  0.79R)  so  that  any  yield- 

6  xd 

depth  combination  for  which  M  =  0.79  should  involve  bubble  behavior 

and  explosion  product  loss  mechanisms  similar  to  those  of  condition 

B.  The  radial-flow  model  fails  for  condition  B  (r^  =  0.99;  implied 

downward  migration)  so  that  the  use  of  equality  of  r^  as  a  sealing 

criterion  is  not  justified  in  this  case.  For  conditions  C,  we  obtain 

M  =  1.^7  and  r_  =  1.20;  both  of  these  criteria  give  virtually  identi- 
6  X 

cal  results  for  yields  greater  than  about  3 ,000  lb.  Results  of  these 

approximate  scaling  calculations  are  shown  in  Fig.  5-24.  For  example, 

7 

if  we  consider  a  nuclear  explosion  of  TNT  equivalent  bubble  yield  10' 
lb.  (5  KT  TNT)  the  relevant  behavior  of  the  bubble  and  the  loss  mech¬ 
anisms  should  be  similar  to  those  of  condition  B  for  a  depth  of  burst 
of  3450 feet, and  to  those  of  condition  C  for  a  burst  depth  of  1300  feet. 

It  should  be  pointed  out  that  the  distributions  of  explosion 
products  found  at  the  second  bubble  maxima  for  these  three  shot  condi¬ 
tions  are  not  necessarily  the  distributions  to  be  expected  at  later 
times.  For  example,  in  the  non-migrating  case,  the  visible  contami¬ 
nated  patch  is  observed  to  grow  after  this  time  by  turbulent  and  par¬ 
ticulate  diffusion,  as  can  be  seen  in  Figs.  5-25  and  5-26.  When  migra¬ 
tion  occurs,  further  vertical  re-distribution  will  take  place  due  to 
"afterflow"  in  the  water  caused  by  the  passage  of  the  bubble.  For 
condition  B  this  can  be  seen  qualitatively  from  the  camera  1  film 
records  in  Fig.  5-27  as  a  vesical  "smearing"  of  the  contaminated 
region  originally  just  below  the  bubble  at  the  second  maximum.  In 
an  attempt  to  qualitatively  evaluate  this  mechanism  for  condition  C, 
four  samples  were  taken  at  intervals  along  the  vertical  explosion 
axis  at  each  of  two  later  times  -  one  shortly  after  the  time  of  the 
third  bubble  maximum  which  occurs  just  below  the  surface  (415  msec) 
and  one  after  the  final  collapse  has  occurred,  and  the  resultant 
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above  surface  "plumes"  are  rising  (63O  msec).  Results  are  presented 
in  Table  4-5  and  Figs.  5-28  and  5-29*  As  can  be  seen,  this  "afterflow" 
apparently  does  cause  substantial  upward  translation  of  the  lost  explo¬ 
sion  products  even  at  these  relatively  early  times. 

CONCLUSIONS 

In  general,  the  most  significant  result  of  ^.lis  study  is  that 
substantial  quantities  of  explosion  products  are  transferred  to  the 
water  environment  during  the  bubble  migration  phase  of  the  underwater 
explosion  event.  Over  the  range  of  migration  strengths  studied 
(which  covers  most  of  the  "very  deep"  ejqilosion  category),  at  least 
3/4  of  this  debris  has  been  lost  from  the  bubble  by  the  middle  of  the 
second  cycle.  For  sufficiently  deep  bursts,  (D  «  60W^^  in  the  kiloton 
range)  as  much  as  95$  of  the  material  may  be  ejected  by  this  time. 

It  is  certainly  not  unreasonable  to  assume  that  additional  material 
will  be  lost  at  subsequent  minima.  Consequently,  it  appears  that  con¬ 
siderable  suppression  of  the  radiological  effects  encountered  at  the 
surface  may  be  obtained  through  these  mechanisms.  The  reservation 
must  always  be  made,  however,  that  material  deposited  in  the  water  may 
eventually  appear  at  the  water  surface  at  later  times  through  diffu¬ 
sion  and  afterflow  of  fluid  generated  by  the  upward  motion  of  the 
bubble. 
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Fig.  5-1  Bubble  Motion  -  Condit: 
(Camera  2  record) 
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107 . 5  msec  -  Second 
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117.5  msec  -  Third 
Maximum 


Bubble  Motion  -  Condition  A 
(Camera  2  record) 
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GOLD  CONCENTRATION  (^o/cm 


SAMPLE  RAOIUS/LOCAL  BUBBLE  RADIUS 


Fig.  5-3  Explosion  Product  Concentration  as  a  Function  of  Sample 
Radius /Bubble  Radius  -  Condition  A. 
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REDUCED  EXPLOSION  PRODUCT  CONCENTRAT 


REDUCLD  RADIUS  (R/o) 


Fig.  5-4  Reduced  Explosion  Product  Concentration  vs.  Reduced 

Radius  at  the  First  and  "Equivalent  Spherical"  Second 
Bubble  Maxima  -  Condition  A 
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Fig.  5 


REDUCED  RADIUS  (R/o) 


5  Fraction  of  the  Total  Device  lying  Within  a  Reduced 
Radius  (R/oc)  of  the  Explosion  Point  as  a  Function  of  Re 
duced  Radius  for  the  First  and  Second  Bubble  Maxima  - 
Condition  A 
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Fig.  5-6  The  Bubble  at  the  Second  Maximum  -  Condition  A 

(Camera  1  record) 
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Fig.  5-7  Transport  Mechanisms  -  Condition  A 
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Fig*  5-3  Bubble  Motion  -  Condition  B 
(Camera  2  record) 
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Fig.  5-3 


Bubble  Motion  -  Condition 
(Camera  2  record) 
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REDUCED  HEIGHT  ABOVE  EXPLOSION  ''TINT  (Z/i) 


Fig.  5-10  Reduced  Explosion  Product  Concentration  Along 
Vertical  Axis  -  Condition  B 
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SAMPLE  RADIUS/CYLINDRICAL  BUBBLE  RADIUS 


Fig.  5-H  Reduced  Concentration  vs.  Reduced  Horizontal 
Radius  -  Condition  B  (5  <  2/a  <  15) 
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Fig.  5-14 
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Fig*  5-15  Bubble  Motion  -  Condition  C 
(Camera  2  Record) 
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Fig.  5-15  Bubble  Motion  -  Condition  C 
(Camera  2  Record) 
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Fig.  5-17 


The  Internal  Density  Structure  of  the  Babble  at 
the  Second  Maximum  -  Condition  C 
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Fig.  5-20  The  Bubble  at  the  Second  Maximum  -  Condition  C 

(Camera  1  record) 
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INCREASING  PRESSURE  (ALL  TO  SAME  SCALE) 
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5-21  Typical  Bubble  Pulse  Wavefoms 
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^6'  5-23  The  Distribution  of  Sq;j 
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Fig*  5-23  The  Distribution  of  Explosion  Products  at  the  First  Maximum 
and  at  the  Second  Maxima  for  Conditions  A,  B,  and  C 
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Fig.  5-24  Shot  Configurations  Characterised  by  Bubble  Behavior  Similar  to  that 
of  Conditions  B  and  C  (Criteria  as  Noted) 
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Fig.  5-26  Equivalent  Spherical  Radii  vs  Time  for  The  Bubble  and 
the  Visible  Batch  of  Explosion  Debris  -  Condition  A 
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(25  msec  After  Third  Minimum) 
Vortex  Ring  has  Formed 


300  msec 

Vortex  Ring  Rising 


Vortex  R: 


Fig-  5-27  The  Effect  of  Afterflow  on  the  Explosion  Product  Distribution  -  Condition  B 

(Camera  1  Record) 
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Fig.  3-28  Late-Time  Effects  -  Condition  C 
(Camera  2  Record) 
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Fig.  5-29  The  Explosion  Product  Concentration 
Along  the  Centerline  Axis  at  Late 
Times  -  Condition  C 
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APPENDIX  A 


THE  SCALING  OF  BUBBLE  BEHAVIOR 

BASIC  CONSIDERATIONS 

In  scaling  the  various  aspects  of  underwater  explosions,  several 
types  of  sealing  are  of  interest.  It  is,  in  practical  terms,  impos¬ 
sible  to  scale  precisely  all  phases  of  the  event  simultaneously.  What 
is  done  in  practice  is  to  confine  attention  to  those  aspects  of  the 
phenomenon  which  are  relevant  to  the  problem  of  interest,  and  then  to 
scale  these  as  closely  as  possible.  For  example,  if  we  were  interested 
in  shock  damage  to  underwater  structures,  we  would  attempt  to  scale 
shockwave  phenomena,  that  is,  to  achieve  Mach  similitude. 

The  general  scaling  requirements  which  are  applicable  to  under¬ 
water  explosions  are  as  follows: 

Similitude  of  Gravitational  Effects 


If  it  is  desired  to  scale  properly  the  effects  of  gravity  on  the 
system,  it  is  first  necessary  to  define  some  ’characteristic  force' 
for  the  model  and  the  prototype.  Then,  similitude  will  be  accomplished 
if  the  Froude  number,  or  rather  the  ratio  of  the  gravitational  force 
to  the  characteristic  force,  is  the  same  for  both  model  and  proto¬ 
type.  That  is. 


If 

jFgrar 

grav 

Fchar 

m 

\ char^ 
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This  nay  readily  be  reduced  to  the  more  familiar  form: 


Where 

v  =  a  characteristic  velocity 
g  -  acceleration  of  gravity 
1  =  a  characteristic  length 
m  denotes  the  model 

p  denotes  the  prototype 


This  requires  that  the  Thomas  number  for  model  and  prototype  be 
the  same,  that  is,  that 

'  P  \  _  /  P  ' 

,  ^vapor  j  m  |  *vapor  j  p 

vhere  P  is  a  characteristic  pressure. 

' imilituae  of  Effects  of  Viscosity  requires  equal  Reynolds  numbers: 

vL\  _  jvL 
v  jm  1  v  P 

vhere  v  is  the  kinematic  viscosity. 

Similitude  of  Cospressibility  (Shockwave)  Effects  requires  equal  Mach 
numbers: 


12:5 


where  c  is  acoustic  velocity 

Other  relationships  of  this  sort  may  he  fabricated  as  needed. 

APPLICATIONS  TO  UNDERWATER  EXPLOSIONS 


In  chapter  2,  th*  *elationships  defining  the  parameters  of  under¬ 
water  explosion  bubble ,  were  presented  assuming  standard  environmental 
conditions,  that  is,  air  pressure  =  1  atmosphere,  the  force  of  gravity 

O 

=  1  the  density  of  the  liquid  environment  =  1  gram/cm,  etc.  In 
the  present  discussion,  greater  generality  is  required,  but  as  will  be 
seen,  the  expressions  used  here  will  reduce  to  more  familiar  forms  if 
standard  conditions  hold. 


Geometrical  Scaling 


As  regards  scaling  of  underwater  explosions,  it  seems  apparent 
that  the  firsn  parameter  which  it  is  desirable  to  scale  is  the  basic 
geometry  characterizing  the  event,  that  is. 


[V 

ni 

m 

hi 

(A.l) 


where  L.  and  L.  are  lengths  characterizing  the  shot.  In  particular, 
^  J 

if  we  take  the  first  maximum  bubble  radius,  or  rather  a  ,  and  the 
shot  depth,  we  can  then  require 


D  .D  (A. 2) 

am-  a  p 


or. 
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where 


D  is  the  shot  depth  in  feet 

a  =  (w/z)1^ 

Z  =  pfg2  +  53  iA 


g 


is  the  fluid  density  in  grams/cm^ 
is  the  acceleration  of  the  system  in  gravities 
is  the  air  pressure  over  the  water  in  atmosphere. 


If  the  sea -bottom  is  near  enough  that  it  will  affect  the  bubble 
motion,  we  must  also  keep  this  distance  geometrically  similar: 


D  -  D 
w 

°  1 

= 

m 

D  \ 


a  j  p 


where  D  =  total  water  depth, 
w 


(A. 3) 


Gravitational  Scaling 

The  driving  force  for  the  bubble  migration  is  the  buoyancy  of 
the  bubble.  Hence,  we  must  reproduce  the  gravitational  effects  upon 
the  bubble  to  preserve  similitude.  Therefore  we  invoke  the  Froude 
criterion: 


IF  ,  \ 

char  •, 

Fchar \ 

F  ! 

grav  m 

,  Fgrav  j 

NoW  Fchar  mchar  ^char/C^char)^ 


and 


L  . 
char 

^char 


a 

T 


(A.4) 
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1 


and,  combining  (9)  and  (10), 


D  ^  Z  /  g  ?  /  p  ?  ,Z  \  ..  Z  p  g  D 

p  p  I °m|  i Kmi  '  p  ,  or  rather  p  p  p  p 

D  3  2  ’l^  ,Z„  Zm  '  »ADm 

in  m  ^ 

we  obtain, 

Pp  ®d  ^p  ®p  +  ^A_  by  definition 

0  g_  D  “  Z  =  p  K  D  +  33  P/rP 
ni  in  m  Hm  “m  m  A 

in 


(A. 11) 


(A. 12) 


Therefore, 

PA  C  g  D  S 

P  _  P  °P  P  _  _» 

P.  p  g  D  Z 

A  x  Tn  m  m 

a 

We  may  now  rewrite  (9)  as: 


(A. 13) 


(A. 14) 


(A.15) 


Thus,  equations  (A.l4)  and  (A. 15),  if  satisfied,  will  satisfy  the 
basic  scaling  requiren-'^.vs  (A.l)  and  (A. 4).  If  we  consider  all  pro¬ 
totype  shots  to  be  fired  in  water,  at  1  atmosphere  pressure,  and  in 
a  gravity  field  of  1  g,  we  may  drop  subscripts  and  wri^e: 

Yield  scale  factor  =  ^  (A.io 

m  ?. 

A 
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L  Dg 

Linear  scale  factor  = 

L  P. 
m  A 


Time  scale  factor  =  -2.  =  2— 


m 


1/2 


(A.17) 


(A.l8) 


Thus,  the  scaling  of  migration  phenomena  may  be  accomplished  by  vary¬ 
ing  the  air  pressure  over  the  water,  the  density  of  the  fluid  and/or 
the  force  of  gravity  in  the  appropriate  way. 


The  Effects  of  Vaporization  and  Condensation 

In  order  to  reproduce  properly  the  effects  of  vaporization  and 
condensation  of  the  environmental  fluid,  it  is  clearly  necessary  to 
reproduce  the  composition  of  the  bubble  atmosphere  -  that  is,  if  a 
point-source  prototype  is  of  interest,  the  model  explosion  must  also 
be  point-source,  so  that  both  model  and  prototype  bubbles  will  contain 
steam. 


If  the  effects  of  condensation  and  vaporization  of  the  bubble 
atmosphere  are  considered  to  be  important,  the  Thomas  criterion  must 
be  adhered  to,  that  is: 


/Pchar  \ 

|Pchar  1 

P 

vapory  n 

P 

\  vapor / 

(A. 19) 


where  Pc^r  is  a  characteristic  pressure  relevant  to  the  phase-change 
process. 


If  we  take  the  pressure  difference  across  the  bubble  interface 

as  our  P  .  ,  we  are  left  with  two  distinctly  different  situations, 

char 

First,  if  we  consider  the  bubble  near  one  of  its  maxima,  we  recall 
that,  near  the  bubble  maximum,  the  internal  pressure  is  a  great  deal 
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less  than  the  hydrostatic  head  Z,  so 

Pchar  =  2  “  Pi  (max)  (A. 20) 

«  Z 

where  P^  is  the  internal  bubble  pressure. 

Near  bubble  minima,  on  the  other  hand,  the  reverse  is  true;  that 
is,  for  reasonable  shot  depths,  (less  than  25,000  feet  or  so),  the 
internal  bubble  pressure  is  a  great  deal  higher  than  the  hydrostatic 
pressure,  so 

Fchar  (miD)  =  Pi  (min)  “  2  (A. 21) 

^  P  (min) 

Considering  first  the  situation  near  maxima,  equation  (A.19) 
and  (A. 20)  require  that 

Zm  ^vapor^  m  ^  PA  +  p8  Dm  (A. 22) 

33  +  dp 

That  is,  the  situation  must  be  arranged  such  that  the  hydrostatic  pres¬ 
sures  at  the  shot  points  for  model  and  prototype  must  be  in  the  same 
proportion  as  the  vapor  pressures  of  the  environmental  fluids.  If 
the  model  explosion  is  fired  in  vater,  this  means,  essentially, 
that: 


Z  =  Z  (A. 23) 

m  p 

if  the  vater  in  the  model  and  prototype  situations  is  at  the  same 
temperature . 

If  ve  now  must  satisfy  our  most  basic  requirements  of  geometrical 
and  gravitational  scaling,  ve  must  then  say: 
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=  1  atmosphere 


v  -  (g i,y  = 


g3. 


since  we  are  now  dealing  with  water. 


If  this  criterion  is  not  satisfied,  the  results  may  he  viewed  as  fol¬ 
lows.  Equation  (A. 22)  implies  that,  if  the  pressure  at  the  shot  point 
is  reduced,  the  prototype  vapor  pressure  is  increased  -  that  is,  near 
the  bubble  maximum,  we  are  scaling  a  large  explosion  fired  in  a  "warm* 
ocean.  In  particular,  if  the  air  pressure  is  reduced  to  the  vapor 
pressure  of  the  water,  we  are  scaling  a  large  explosion  in  a  boiling 
ocean,  and  the  pressure  may  be  reduced  no  further.  Physically,  the 
significance  of  this  is  the  observed  boiling  of  the  babble  interface 

near  the  bubble  maximum  for  shots  fired  at  hydrostatic  pressures  near 

9  19 

the  vapor  pressure  of  water.  The  effect  of  such  boiling  is  not 

entirely  clear. 


The  problem  of  reproducing  vaporization  and  condensation  phenomena 
near  the  bubble  minimum  is  somewhat  different.  As  we  have  seen  (Sq. 
A.21). 


Pchar  (min)  =  Pi  (min) 

where  P^  (min)  is  the  peak  internal  pressure  attained  by  the  bubble 
near  its  minimum  size.  The  problem  then  becomes  that  of  evaluating 
P^  (min).  Ref.  9  presents  the  experimental  result  that  for  bubbles 
which  do  not  migrate  appreciably,  the  volume  of  the  bubble  minimum 
is  approximately  ten  times  that  of  the  "initial  bubble"  (or,  the  "equi¬ 
valent  charge  volume"),  and  that  this  volume,  in  contrast  to  the  maxi¬ 
mum  bubble  volume,  is  apparently  independent  of  the  loctl  hydrostatic 
pressure . 

Consequently,  if  the  bubble  atmospheres  of  model  and  prototype 
are  composed  of  similar  materials,  such  that  the  ratio  of  specific 
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heats  (7)  is  the  same  for  both,  then,  since  the  expansion  and  contrac¬ 
tion  are  adiabatic, 

3?i(0)  Vj  =  Pi(min)  VmJ  (A. 2k) 


vhere, 

P^Co)  =  pressure  in  initial  bubble 

V  -  volume  of  initial  bubble 
o 

V  .  =  volume  of  bubble  minimum 

min 

Hence,  for  the  nonmigrating  case, 

PjL(min)  «  Pi(0)/107  (A. 25) 

Nov,  since 

R±  «  0.135  Wl/3  feet  (A. 26) 

and  the  total  internal  energy  in  the  initial  bubble  is  just  the  bubble 
energy, 

P±(0)  s  constant, 

so  that  the  internal  pressure  at  the  bubble  minimum  is  independent  of 
both  yield  and  hydrostatic  pressure  for  the  non-migrating  case. 

For  the  migrating  bubbles,  the  volume  of  gas  in  the  bubble  minimum 
is  not  veil  defined,  since  in  all  likelihood  most  of  the  "visible" 

6  9 

volume  of  the  minimum  is  occupied  by  the  water  in  the  undervater  jet.  ’ 
It  is  certainly  not  unreasonable  to  assume  that: 

*KbU)  <»s>  -  Ei(mln)  <°>  *  f '  <V 

vhere  f 1  is  an  unspecified  function  of  the  specific  migration,  which 

21 

is  equal  to  1  for  non-migrating  bubbles.  If  this  is  the  case,  t^en 
the  internal  pressure  at  the  minimum,  P^(min)  is  e  function  only  of 
the  strength  of  migration.  Consequently,  the  Thomas  criterion  will 
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be  satisfied,  and  evaporative  and  condensive  phenomena  will  be  scaled 
near  the  bubble  minimum,  if  the  requirements  of  gravitational  and 
geometric  scaling  (which  define  the  migration)  are  satisfied. 


(A.  31) 


Multiplying  equation  (A. 3})  by  1  on  either  side  (that  is,  by  equation 
(A. 30),  we  obtain. 


$.  -  p  i/2 
W  "  A 


(A.  32) 


For  example,  if  W 


10  W  ,  then  this  means 
m 


=  100  atmospheres 

Thus,  to  satisfy  equation  (A. 9)  for  gravitational  scaling, 

=  1000  gravities 

That  is,  ve  are  actually  requiring  that  the  air  pressuic  over  the 
water  be  increased,  leaving  only  g  as  an  independent  variable.  Clearly, 
Reynold's  similitude  is  not  feasible. 


It  should  be  pointed  out,  however,  that  viscosity  is  not  important 
for  Reynold's  numbers  greater  than  about  100*  or  so  -  we  then  enter 
the  regime  of  turbulent  drag  forces.  Reynold's  numbers  are  always 
very  large  during  the  bubble  migration  phase  of  an  underwater  explo¬ 
sion.  The  question  is,  then,  how  are  turbulent  drag  forces  reproduced? 
As  was  done  previously,  we  set: 


I  turb1 
Pchar  I  m 


(F 


turb 
char  /  p 


Now,  the  drag  force  may  be  characterized  as  follows: 

2 


turb 


CD  *  A  X  pf  V 


where 


A  =  cross  sectional  area 
^f  =  density  of  medium 
v  =  velocity  of  body 


(A.33) 


(A- 3*0 
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(A. 35) 


=  drag  coefficient 

So,  a  characteristic  drag  force  is  as  follows: 


Fturb 


Pf  £ 
2  2 


And,  as  before,  our  general  characteristic  force  is: 

4 

F  =  pf  51- 

char  i  2 


(A. 36) 


Consequently,  if  the  configuration  of  the  bubble  and  the  bubble  migra¬ 
tion  is  reproduced,  turbulent  drag  forces  will  scale. 

The  Effects  of  Compressibility  -  Shock  Effects 

In  order  to  scale  the  effects  of  the  shockwave,  we  must  set  the 
Mach  number  for  model  and  prototype  equal,  that  is. 


Vchai|  |Vcharl 

(A.  37) 

c  Jn  "  i  c  ip 

where 


c  =  acoustic  velocity  in  the  medium 

_  71/2 


vchar  =  •  2 


So  the  requirement  becomes: 


If  this  criterion  is  to  be  satisfied,  the  air  pressure  in  the  model 
environment  must  remain  at  one  atmosphere. 


As  we  have  seen,  however,  the  interaction  between  the  shockwave 
and  the  bubble  ends  at  very  early  times  in  the  first  cycle,  before 
any  migration  has  occurred.  Hence,  for  purposes  of  studying  bubble 
behavior,  Mach  similitude  is  not  required. 
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APPENDIX  B 


LIST  OF  SYMBOLS 


A  =  Cross-sectional  area 

c  =  Acoustic  velocity 

C  =  Explosion  product  concentration 

C__  =  Background  explosion  product  concentration 

iJtr 

Cp  =  Drag  coefficient 

2=0,=  Shot  depth 

D^  =  Depth  of  second  bubble  maximum 

D^  =  Depth  of  vauer 

Eg  =  Shockwave  energy 

Ep  =  Total  energy 

g  =  Acceleration  of  gravity 

k  =  W/W* 

L  =  Lchar  =  characteristic  length 
Mg  =  Specific  migration  =  (D^-D^/R^ 

Subscript  m  -  Applies  to  model 
Subscript  p  -  Applies  to  prototype 
p  =  Momentum 

F  =  ^char  =  Characteristic  pressure 

PA  =  Air  pressure  in  atmospheres 

P^  =  Internal  bubble  pressure 

Pygpor  =  Vapor  pressure 

r  =  Radius  of  sample  from  ejqplosion  axis. 

rT  =  Period  ratio  =  \f2^  =  0)/^^  --  180') 

R  =  R^  =  First  maximum  bubble  radius 
Ri  =  Radius  of  initial  bubble 
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R  =  R(«/2)  =  First  maximum  side  bubble  radius 
s 

R(t)  =  Bubble  radius  as  a  function  of  time 
R^  -  Second  cycle  maximum  bubble  radius 

S  =  S(p)  =  Surface  correction  coefficient  (a  function  of  p) 

T  =  First  bubble  cycle  period 

^l/2  =  cyds  "half  period"  (time  to  maximum  expansion) 

T  =  T(it/2)  =  First  side  bubble  cycle  period 
s 

v  =  Characteristic  velocity 

W  =  TNT  equivalent  bubble  yield 

W*  =  Charge  weight  or  total  released  energy. 

2  =  Altitude  of  sample  above  explosion  point 
Z  =  Z^  =  Hydrostatic  pressure  at  shot  point  in  feet  of  water 
a  =  Characteristic  shot  length  =  (w/z)^3 
P  =  Geometrically  scaled  depth  =  D/a 
y  =  Ratio  of  specific  heat  (C^/Cy) 


AD12  =  D,  -  D2 

e  =  *he 

pd/o3  =  Reduced  unit  of  explosion  product  concentration  -  "micro-devices 
per  cubic  alpha" 

V  =  Kinematic  viscosity 

pe  =  Mass  density  of  explosive 

p  =  pf  =  Density  of  fluid 

T  =  Characteristic  shot  time  =  W^3p^^/Z*^ 

0  =  Bubble  latitude 

cp  =  Energy  density  of  explosive 
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